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IV. 1 PROPULSION 


Propulsion systam parformanca is computad during angina sizing and 
whenever aircraft parformanca is computed. Iwo separata sets of propulsion 
subroutines are used for jet eind propeller driven aircraft as illustrated 
in Figtira IV. 1.1. 

TV. 1.1 Turbojet/Fan Propulsion Subroutines 
The turbofan/turbojet engine performance methodology is based on tabu- 
lated performance data for specific engine cycles. Currently, seven different 
engine cycles are available representing a wide range of operational and 
conceptual engines. Data for these engines are contained in subroutines 
ENGDTl-7. Tabular engine data for turbofan/ turbo jet engines may also be input 
in which case ENGDTT is used to determine engine performance. 

Performance data for each of the engine cycles may be scaled up or down 
to simulate an engine of arbitrary size. Engine size is expressed in terms 
of sea level static airflow zuid is determined in subroutine ENGSZ. Engine 
performance is determined by subroutine ENGINE as a function of the flight 
altitude, Mach number, and engine power setting. 

ENGSZ determines the engine size necessary to meet selected performance 
requirements. The engine size is expressed in terms of rated /flow under 
normal sea level static conditions. The engine is first sizv?.' match 
cruise drag, with an option to specify a rate of climb margin. An option is 
then provided to resize the engine so as to match a required takeoff distance 
or to match one-ongine-out requirements on the aircraft rate of climb (PAR, 
Part 23 or Pzort 25.). 
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T-AI 


TUPBOFAN 

PROPELLEi; 

ENGSZ (510) DETERMINES ENGINE SIZE 

ENGSZ (615) DETERMINES r^lZB 

ENGINE (110) ENGINE PERFORMANCE 

ENGINE (480) POWER PLANT PERFORMANCE 

ENOiTl (50) ENGINE DATA 

PWRPLT (65) RECIPROCATING ENGINE PERFORMANCE 

ENGDT2 (80) ENGINE DATA 

TURBEG (2”0) TURBOPROP ENGINE PERFORMANCE 

ENGDT3 (70) ENGINE DATA 

ENO)AT (120) PROPELLER CHARACTERISTICS 

ENGDT4 (70) ENGINE DATA 

PERFH (490) PROPELLER PERFORMANCE 

ENGDT5 (70) ENGINE DATA 

COST (35) PROPELLER COST 

ENGDT6 (70) ENGINE DATA 

GEARBX (45) GEARBOX COST, WEIGHT 

ENGDT7 (145) ENGINE DATA 

VZAIT (30) PROPELIER WEIGHT 

ENGDTT (35) ENGINE DATA 


NACDG (35) NACELLE DRAG 

PNOYS (50) CCWTROLS PROP/ENGINE NOISE CALCULATION 


ZNENG (55) ENGINE NOISE CHARACTERISTICS 


ZNOISE (105) PROPELLER NOISE 


FIGURE IV. 1.1 PROPULSICW SUBROUTINES 



Tha angina aizing problam is itarativa/ because of the affect of nacalla 
gaoraatry on total aircraft drag. The input flag KNAC accounts for nacelle 
drag by tha following means: 

KNAC « 0 Nacelle drag is included as an angina performance 

penalty. Nacalla size is a function of angina size 

1 

and is computed during engine sizing. 

KNAC »« 1 Nacelle drag is accounted for as an aerodynamic force. 

Nacelle size is a function of engine size and is 
computed during engine sizing 

KNAC « 2 Same as KNAC » 1 , but nacelle dimensions are input 
and remain fixed. 

When KNAC ■ 0 or 1/ the nacelle size is initially estimated as a function of 
aircraft gross weight since the engine size is uiiknown. Tha required engine 
size is then computed, and, based on this engine size, an improved estimate of 
the nacelle dimensions is made. The nacelle diameter is computed from the 
sea level statis airflow (WASLS) , the feui face Mach number <SM1D, input) , and 
the fan hub-to-tip ratio (HBTP, input) using one-dimensional isentropic 
compressible flow theory; the nacelle length is computed from the diameter 
and cin input nacelle fineness ratio (XLQDE) . Based on these dimensions the 
nacelle drag is recomputed, and the engine sizing process is repeated once. 

When KNAC * 0, the nacelle drag is computed exactly as when KNAC » 1. 
Engine specific thrust 2 ind specific fuel consumption are adjusted for nacelle 
drag in subroutine NACDG. 

The flight condition input flag for engine sizing (JENGSZ) can take on 
the following values; 

^ Applies to tirrbof an/turbo jet configuration 
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JBNGSZ 


FLIGHT CONDITION 


0 Sis* »ngin« for cruisa flight condition 

1 Size engine for cruisa and taksoff flight conditions 

2 Size angina for cruisa and takaoff and ona-angina- 

out climb flight conditions 

3 Siza angina for cruisa and ona-angina-out climb flight 
conditions 

4 Engina thrust is spacifiad. 

Tha anginas ara initially sized at tha design cruisa flight conditions except 
whan tha angina siza is input (JENGSZ » 4) . This means that/ at cmisa 
power setting/ tha anginas must produce total thrust equal to tha cruise drag 
of the aircraft. If a cruisa climb margin is specified (ECCRU) / the engines 
must also have enough excess cruisa thrust to meet this margin. 

The required engine size, expressed as sea level static airflow (WASLS) 
is ccanputed by scaling tha performance of the reference engine to match tha 
required cruise tiirust. Engine performance is scaled by assuming that at a 
given altitude, Mach number, and engine power setting tha specific thrust 
(SPN * thrust per unit airflow) and percent corrected airflow (PCWAC « corrected 
airflow/A^ASLS) of the scaled engine are the same as for the reference engine. 
Thus, the sea level static airflow of the engine is computed by 
Cruisa Airflow » required thrust/SFN 

Corrected Cruise Airflow >■ 

^ruisa Airflow x 

WASLS » Corrected Cruise Airflow/PCWAC 


'' Total tWp 

SLS Temp 

Total Pressure 
SLS Pressure 
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Wh«n JKNGSZ « 1 or 2, th« taka-off diatanca of tha aircraft (with angina* 


sized for cruise) is computed (ENGSZ call* PERPRM which call* TAKOPF) and 
compared with the input required distance (XTORQ) , tlie required t^lke-off 
distance may be for high altitude and for hot day conditions. If the computed 
distance exceeds tha required distance, then the engines are resized by 
adjusting the airflow to meet this requirement. 

Federal Air Flegulations Parts 23 and 25 establish climb requirements. 

For example, FAR Part 25 requirements eire surranarized in Figure IV. 1.2. When 
JENGSZ is input as 2 or 3, ENGSZ computes the climb performance of the air- 
craft in accordance with Part 23 or 25 and compares the computed performance 
with the required performance. If necessary, the ,jines are resized so that 
the aircraft meets the most critical requ'rement. 

ENGSZ includes an option for sizing the engines for an input turning 
performance requirement. This option (JTRSZ « 1) must be used in conjunction 
with one of the engine sizing options described above (JENGSZ « 0-3; may not 
be used with JENGSZ » 4) . 

The user must specify the required turn load factor, altitude and Mach 
number (XLFTRN, HTURN, EMTURN - input) . ENGSZ computes the thrust required 
to execute the turn and compares this thrust with that available from the 
engines at the desired power setting (engines as sized for cruise, takeoff, 
or climb) . If insufficient thrust is available, engine thrust is set equal 
to that required, ^md a new sea level static airflow is determined. An 
additional iteration is performed to account for resized nacelles when nacelle 
size is a function of engine airflow (KNAG •* 0, 1). 

Turning performance may be limited by the maiximum lift coefficient in 
tha turn configuration. This limiting lift coefficient may be specified by 
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the vLsar (input as CLTLMT; default » 1.0). If tlxa required lift coefficient 
in the tvim exceeds the maximua turn lift coefficient/ the turn load factor 
is automatically reduced to the value achievable by the aircraft at itr 
limiting turn lift coefficient. 

The simplest engine sizing option (JENGSZ » 4) is for the user to specify 
the rated sea-level static thrust (THIN/ lbs.) of one engine. In this case 
engine sizing at cruise is bypassed. Engine-out climb performance is canputed 
as when sizing for climb? however, if a climb deficiency is detected, the 
engines are not resized. 

When the engine size is input, several additional inputs are required. 

The nacelle size must be specified (KNAC « 2, EIjN and DBARN) . Ixa addition, 
the engine, nacelle, and pylon weights (WENG, WNAC, WPYLON) must be input if 
non-zero values are desired. 

IV. 1.1.1 Engine Performance at a Specified Flight Condition - Subroutine 
ENGINE . Subroutine ENGINE determines engine performance at a specified 
altitude and flight Mach number. Engine performance is described by thrust, 
airflow, fuel flow, specific thrust, per cent corrected airflow, and thrust 
specific fuel consumption. Performance data for different engine cycles are 
contained in subroutines ENGDTl-7, eind ENGDTT as functions of altitude, Mach 
number, and power setting. Either power setting or required thrust is 
specified, and the other is to be found. 

The perfoirmance of a particulair engine may be scaled up or down by 
assuming that at a given Mach number, altitude, and engine power setting, the 
specific thrust (SEN » lbs. thrust/lbs. /sec airflow) , the specific fuel 
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CONFIGURATION 

IRQ 

ALT. ABOVE 
AIRPORT, FT 

TAKEOFF FLAPS. 
LANDING GEAR DOWN 
ONE ENGINE OUT 
(FIRST SEGMENT) 

1 

0 

TAKEOFF FLAPS 
LANDING GEAR UP 
ONE ENGINE OUT 
(SECOND SEGMENT 

2 

250 

CLEAN CONFIGU- 
RATION 

ONE ENGINE OUT 
(FINAL TAKEOFF) 

3 

1500 

APPROACH CONFIGU- 
RATION 

ONE EliGINE OUT 

4 

0 

LANDING CONFIGU- 
RATICXt 
ALL ENGINES 



5 
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FIGURE IV. 1.2 


ENGSZ; FIRST 


CLIMB VELOCITY 
(Vc> (KNOTS) 


V » V 
C LOF 

V < 1.2 V 

LOF" STALL 


^ ^STALL 


Vc<1.5*Vstall ® 

Vc4l.l*VsTALL ® LANE 


* "Wall ® ^ 


REQUIRED CLIMB 
GRADIENT 

SOURCE 

POSITIVE 
(Rt 1 FT/S) 

FAR 

25.121(a) 
2 engine 

; 1 

1 2 . 4 % 

^ Vci^xmb*‘024 
F/S) 

FAR 

25.121(b) 
2 engine 

1.2% ^ 

(RC» ^clIMB 

F/S) 

FAR 1 

25.121(c) 1 
2 engine : 

2.1% 

(RC ^ VclIMB**®21 
F/S) 

FAR 

25.121(d) : 
2 engine 

3.2% 

Vclihb*-032 

F/S) 

FAR 
25.119 
i 2 engine 


D SECOND SEGMENT CLIMB REQUIREMENTS (FAR PART 25) 
































conautaption (SFC - Iba. of fu«l par hour par pound of thruat) , and tha par 
cant corractad airflow (correctad airflow dividad by aaa laval atatic air- 
flow) of tha acalad angina ara the aans aa for tha unacalad angina. Once tha 
aea level atatic airflow (WASLS) la eatabliahed by subroutine ENGSZ, ti'-a 
scaled engine performance at the aneclfied operating point follows immediatelyj 

Airflow: Wg “ WASLS x PCWAC X 5/^5* 

Thrust Ffj SFN X Wq 

Fuel Flow: Wp - SFC/F^ 

where 5 ■ total praaaura/SLS praaaura 

Q » total temperature/SLS tasnparature 
The different functions of ENGINE are con "'^llod by the indicator KENG: 

KENG » 1 (Variable correctad rotor speed) 

This option is used when the required thrust is known, for axaaapla during 
cruise at a specified altitude and Mach number. The engine power setting is 
found at which engine thrust is equal to required thrust. The fuel flow at 
this power setting is used in the range calculation. Engine data art- not 
scaled. 

KENG » 2 (Idle power setting) 

Engine perfomuance at idle la used during tha taxi segment and during the 
landing calculation (only if the engines have already been aired) 

KENG » 3 (Maximum cruise corrected rotor speed) 

This option is used during engine sizing at cruise and during performance 
calculations at maximum cruise power. Engine power setting is known, and 
engine performajjce is determined. If the engines have already been sized, 
engine performance is computed from WASLS, PCWAC, SFN, and SFC. During engine 

IV- 1 
8 



sixlng/ spaclflc aiiigln* parfomanc* (PCWAC, S7C) is uaad tha raquirad cruise 
thrust to computa engine si*e (see Section IV, 1). 

KENG - 4 (Variable Corrected Rotor Speed 

This option is the same as KENG - 1, except that engine data are scaled using 
WASLS, PCWAC, SFN, and SFC. 

KENG « 5 (Maximum Corrected Rotor Speed) 

ENGINE determines the scaled engine perfomumce available at the operating 
flight condition and maximum engine power setting. This option is used during 
take-off and climb. 

KENG « 6 and 7 (Maximum Continuous Power Setting) 

KENG - 7 (Meiximum Climb Power Setting) 

KENG is automatically set according to the type of performance calculation 
whenever ENGINE is called. Take-off, climb, acceleration, and turn performance 
are normally computed at maximum corrected rotor speed with KENG * 5. Cruise 
performance at normal rated power is computed at maximum cruise power (KENG « 
3). Cruise performance at a specified Mach number is computed with KENG ^ 4. 

Takeoff, climb, acceleration or turn performance will be computed at 
maximum continuous (KENG « 6) or maximum climb (KENG •* 7) power settings rather 
than maximum power if the variables KODETO, KODECL, KODEAC, and KODETR are 
input as 6 or 7 (default values for these variables are 5) . 

For all of these options, it is assumed that the following normalized 
engine parameters are the same for the scaled engine as for the model engine i 

1. Percent corrected airflow, PCWAC 

2. Specific tlirust, SFN 

3. Specific fuel consumption, SFC 
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Th« two additional oparational paramatar* naadad ara tha ratio* of total 
prasBura and total tamparatura to sea laval static valuas of thasa parameters, 
in terms of which the "corrected" values of thrust, airflow, and fuel flow 
can be found. 

other input parameters included in tha subroutine call statment are 
FN ■ engine thrust, lb (if power setting is being computed) 

SFC *• engine specific fuel consumption, lb per sec per lb 
FAR " fuel-to-air ratio 

PO «• ambient static pressure, lb per sq ft 
TO * ambient static temperature, deg R 
SMN - engine Mach number 
HN »■ altitude, ft 

Primary output quantities of tha routine ara 

FN - engine thrust, lb (if power setting is knc«m) 

WF - fuel flow, lb per hr 

IV. 1.1. 2 Engine Data Subroutines. Turbofan/ turbojet angina performance 

data is available in the following eight subroutines: 


SUBROUTINE 

ENGINE 

lENGSC 

ENGDTl 

GE CJ610-6 (Turbojet) 

1 

ENGDT2 

Garrett TFE 731-2 (Turbofan) 

2 

ENGDT3 

UACL JT15D-1 (Turbofan) 

3 

ENGDT4 

Lycoming ALF-502 (Turbofan) 

4 

ENGDT5 

GE CF 34 (Turbofan) 

5 

ENGDT6 

GE TF 34 (Turbofan) 

6 

ENGDT7 

Conceptual GE T700/F1 QCGAT 

7 

ENGDTT 

Tabular input engine data 

•«0 
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Th# uppropriftta angina cycla ia aalaotad by tha valua of lENGSC wilaich 1» 
read from the data deck title card (cola. 75-76) . 

Each engine data si:ibroutine contains tabulated parfonaanca data for a 
specific engine cycle. For ENGDTl-7 these data consist of corrected thrust/ 
corrected airflow, and corrected fuel flow tadjulated as functions of engine 
power sotting and flight Mach number. These relationships eire illustrated 
schematically in Figure IV. 1.3. The effect of altitude is contained impli- 
citly in the total temperature and total pressure ratios (5a3id5). 

These normalized or specific quantities, used to scale engine performance, 
are computed each time ENGDT is called: per cent corrected airflow (PCWAC) , 

specific thrust (SFN) , and specific fuel consumption (SFC) . 

Engine power setting is established according to the valua of KENG, which 
is set in subroutine ENGINE: 

KENG >» 5 maximum power 

1, 4 variable ixwer 
3 maximum cruise power 

2 idle power 

6 maximum continuous power 

7 maximum climb power 

Typically, engine power setting is expressed as either the ratio 
of turbine inlet temperature to engine face total temperature (T5/T2) or 
per cent corrected rotor speed (N/V^T/Nmax) • Each engine cycla does not neces- 
sarily possess the complete range of power settings indicated by KENG. For 
example, maximum continuous power may be identical to mar.imum cruise power. 


etc 






ENGDTT dlff«r« from WiGCfTl-7 in that tha angina ceorractad parforraance data ara 
tabulated as explicit functions of altitude as wall as power setting and 
Kach number. In addition, these tables must be rend from cards at the begin- 
ning of each run {see input deck description) . These engine data are read, 
stored, and interpolated with the aid of several special utility subroutines 
and functions {MAPS, ST0RE3, TTABX, TABX, DTABX, BISC) . 

Engine idle performance {KENG • 2) is determined slightly differently 
for each engine cycle. In the simplest cases, numerical values are assumed 
for idle corrected thrust, corrected fuel flow, and corrected airflow. Two 
engine cycles (ENGDT2, 3) express corrected engine-idle performance as expli- 
cit functions of altitude. 

All of these subroutines have the same list of 13 arguments, which 
effectively specify all aspects of engine performance. These arguments are 
FN • engine thrust, lb 

WF » fuel flow, lb per hr 

WA • airflow, lb per hr 

PCWAC ■ percent corrected airflow 

■ corrected airflow/SLS airflow 
SFC » pecific fuel consumption, lb per hr par lb thrust 

SFN » specific thrust, lb per lb per hr airflow 

FAR “ fuel/air ratio, lb per sec of fuel per lb par hr of air 

P2 * static pressure at inlet, lb per sq ft 

T2 ■ static temperature at inlet, deg R 

HN « altitude, ft 
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SM I'light Mach watabar 

PTHROT »• par cent maximum throttle setting 

KENG "• angina power setting indicator (0-7) 

The output quantities are the first seven of these parameters (input being 
the last six), and they are determined by interpolation in the numerical 
performance tables which mabe up the major part of the ENGDATl-7 subroutines. 

IV. 1.1. 3 Nacelle IjOases Computation, Subroutine NACDG . This subroutine 
corrects the engine specific fuel consumption and specific thrust to include 
losses due to nacelle drag. It is called by ENGINE when KNAC « 0, and the 
subroutine subtracts nacelle drag from gross engine thrust to provide a net 
engine thrust. The drag coefficient is approximated by using flat-plate 
turbulent boundary layer theory, and this depends on the cylindrical nacelle 
dimensions and the Reynolds number of the nacelle. 

The input argiunents of the subroutine relate to the flight conditions 
auvl the aircraft geometry, as needed for drag coefficient computation. These 
include static pressure PO, Mach number SMN and wing area SWING. The output 
quantities are the corrected specific fuel consumption and specific thrust 

(SFC, SFN) and the total airflow, w. . 

A 

IV. 1.2 Propeller Propulsion Subroutines 
The propeller propulsion subroutines in GASP are used to simulate the 
performance of reciprocating, rotary combustion, and turboprop propulsion 
systems. Several of these subroutines replace equivalent turbofan engine 
subroutines and thus have the same names . 

Svibroutine ENGS2, like its turbofan equivalent, controls several angina 
sizing options. It determines the engine size necessary to allow the aircraft 
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to »*«t a a«t of parfomaOiC* raquirsaants. Th* emgina parformanca aubroutina, 
ENGINE, i* callad during angina aising and aach tlma propulaion ayatam 
performance ia required.lt relates the performance of the propeller, controlled 
by subroutine ENGDAT, to the perfornuinca of the powerplant, computed in 
aubroutina PWKPLT for rotary cocnbusion and reciprocating anginea and in TURBEG 
for turboprop engines. 

Reciprocating and rotary combustion engine performance (aubroutina 
PWRPLT) is based on generalised, non-dimensional relationships between power, 
engine speed, and altitude. Both normally aspirated and turbo-charged engines 
may be simulated. 

Turboprop engine performance ia based on tabulated data for a specific 
engine cycle. Currently, data for the Garrett TPE 331-1 turboprop are 
included in the program (subroutine TURBEG) . The performance of the baseline 
engine may be scaled up or down. 

Propeller performance is computed as a function of the propeller geometry 
and operating condition using generalized performance relationships. These 
data are contained in subroutine PERFRM, which is called by ENGDAT whenever 
propeller performance is to be computed. 

IV. 1.2.1 Propeller Subroutine, ENGSZ . The propeller subroutine ENGS.^ 
determines the engine size necessary to meet cruise, and, optionally, take-off 
and/or climb requirements. Engine size is expressed as maximum sea level 
horsepower for reciprocating, turboprop, and rotary combustion engines. 

The engine size may be input directly as HPMSLS by also inputting KODECR 
7. In this case, the engine sizing process is bypassed. For fixed pitch 
propeller configurations, propeller blade angle (BLANG) must also bo input 



wh«n tha aiigJLiia «isa is input. Tha input flag JENGSZ datarm.in.as th,a sising 
options, as follows! 

JENGSZ ■ 0: sisa at cruisa only 

1: sisa at cruise and takeoff 

2: size at cruisa, takeoff and climb 

3j size at cruise and climb 

In all casas, tha engines are sized to provide the required cruise thrust 
at an input flight condition and engine operating point. In tha last three 
cases, the engine size may be increased so that the takeoff aid climb perfor- 
mance requirements are met. 

The cruise condition is specified by altitude, Mach number, cruisa 
weight, required cruise rate of climb and engine operating point. When all 
of these parameters axe specified, the thrust required for a given rate of 
climb is easily expressed as a function of cruise drag, weight, velocity, and 
rate of climb. Engine sizing at cruise involves ccmiputing the horsepower 
necessary to produce the required thrust at the design cruise condition. 

Since this cruise power is a specified fraction of maximum sea level power, 
the rated power of the engine may be computed. 

As was discussed in Section IV. 1.1 there exist thr- nacelle drag 
accoiuiting options j however, only two are available for propeller power 
configurations. 

KNAC “ 1 Nacelle drag is computed in ENGSZ as engine size is 

determined. Nacelle size is a function of engine power. 



KHAC ** 2 Nacslla dr*jg lx ooc^putflkd from input nacalla dlmanxion* 

in AERO. Thix option should ba uaad to zero out nacelle drag 
for single engine nose inounted anginas or other configurations 
with buried engines. 

In any case/ nacelle drag is found as an explicit function of nacalla 
Reynolds number and wetted area, euid the watted area and Reynolds number are 
both numerical functions of aircraft geometry. 

The engines may be resized to enable the aircraft to meat an input take- 
off field length requirement (JENGSZ » 1 or 2) . This requirement may include 
high altitude and/or hot day conditions. The take-off performance of the 
aircraft with engines sized for cruise is computed (calls to subroutines PERI’RM) 
and TAKEOFF) and if the computed take-off distance exceeds the required 
distance, ENGSZ iterates on engine power tmtil the take-off requirement is met. 

When JENGSZ » 2 or 3, the program ccxnpares the aircraft's climb performance 
with the requirements established in Federal Aviation Requirements, Part 23 or 
25 (shown in Figure IV. 1.3). If one or more of these requirements is not met, 
the engine size is increased to satisfy the most critical requirement. 

Propeller diameter is an input variable; however, if the engines are 
resized for take-off and/or climb, propeller diameter may or may not be 
changed according to the input variable JSIZE: 

JSIZE « 1 Increase power but leave propeller dizimatar constant 

** 2 Increase both power and diameter but keep propeller disk 
loading (power/area) constant 

When the engines are sized such that all take-off zind climb requirements 
are mat, additional engine resizing may be needed if KNAG "1. In this case 
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th« origins raslslng is rapastad, b«.sad on tha naw astimata of nacalla size. 

And, if angina is resized only at climb, with KNAC *■ 1, tha climb performance 
must be recomputed using the new estimate of nacelle drag. Again, in this 
case, engine power may be increased tmtil desired performance is attained. 

Engine sizing for fixed pitch propeller configurations is handled somewhat 
differently th^ul for constant speed propeller configurations. The engine is 
initially sized at cruise just as for constant spaed configvurations. The 
propeller blade angle computed at cruise is held fixed for subsequent perfor- 
mance calculations. Thus, initial sizing at cruise establishes both the 
engine size ^md the propeller blade angle. 

If the climb performance of a fixed pitch piopallar aircraft with the 
blade emgle set for cruise does not meet all tha climb requirements (JENGSZ *• 

2 or 3), then the blade emgle is decreased by two degree increments, and the 
climb performance is recomputed. The largest blade angle for which all require- 
ments are met is fixed as the new blade 2 mgle. If blade emgle reductions alone 
fail to sufficiently improve climb performance, then engine power is increased 
^uld the propeller blade emgle is set for tha critical climb requirement. 

IV. 1.2. 2 Power Plant/Propeller Matching, Subroutine ENGINE . The most 
important function of subroutine ENGINE is to relate the performance of the 
powerplant (piston engine, rotary ccanbustion engine, or gas turbine engine) to 
the performance of the propeller (or other propulsor, such as a Q-fan) . ENGINE 
is called dviring engine sizing and during performance calculations when either 
propeller thrust or engine power is known, and the other is to be determined. 

Engine deals with several types of propulsion performance problems. The 
flag KENG which is passed to ENGINE through its argument list is used to 
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»p»cify th« typ« of computation r*quir*d. Th* indicator KODE is ««t in 
ENGINE according to the value of KENG and further specifiea the type of 
computation. The types of engine performance calculations are suamarired 
below; 

KENG » 3 (KODE » 1, 2, 3, or 4) 

This option is used during engine sising at cruise whan the 
required propulsive thrust is known, emd the required angina sise 
is computed. Normally, the engine cruise operating point (percent 
power and per cent RPM) is specified (KODE » 4) . Options exist for 
adjusting propeller diameter (KODE * 1 or 2) or the propeller cruise 
RPM (KODE =» 3) to maximize Propeller efficiency at the design cruise 
flight condition. The value of KODE may be specified by inputting 
the appropriate value for KODECR. 

KENG « 1 or 4 (KODE » 5 or 6) 

This option is used during cruise performance calculations when 
the cruise Mach number ^md altitude nxe specified. The engine and 
propeller characteristics are fixed, and the required propulsive 
thrust is known. ENGINE fiivis the power setting required to drive 
the propeller and the corresponding fuel consumption. Propeller 
RPM is either specified at the design cruise value (KODE *• 6) , or 
ENGINE will select the RPM which maximizes propeller efficiency 
(KODE “ 5) . The value of KODE may be selected by inputting the 
appropriate value for KODETH. 

KENG » 0 or 5 (KODE « 7) 

This option is used when the engine sire, power setting, RPM, and 
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aircraft flight spaad ara kr. m and th* raaultant thrust and fual 
consumption are computed. This situation corresponds to aircraft 
equipped with var iaible pitch propellers during tahe-off, climb, 
acceleration, and cruise at a specified power setting. 

(KODE « 8) 

This option is used to predict the performance of aircraft equipped 
with fixed pitch propellers during full throttle operation at a 
specified airspeed (take-off, climb, acceleration) . ENGINE finds 
the engine and propeller RPM at which full throttle power available 
ec[uals the power absorbed by the propeller. Having found the equi- 
librium RPM, ENGINE finds the corresponding propeller thrust and 
engine power and fuel consumption. 

ENGINE finds propeller performance by calling subroutine ENGDAT. Power- 
pl 2 mt performance is found frcan subroutine PWRPLT for piston and roteury 
combustion engines and TURBEG for turboprop engines. 

One important function of subroutine ENGINE is to insure that the 
operating conditions of the propeller and power plauit are compatible. Specifically 
the power required to turn a propeller at some RPM and flight condition must 
not exceed the maximum power available from the engine at that RPM and flight 
condition. If power required exceeds power availeible at the specified RPM, 

ENGINE seeks some other engine speed where power available is sufficient to 
drive the propeller. 
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17.1.2.3 Povwr and gual lFla» Computatlong , gubroutina PWKPLT . Tlii« 


•labroutine ccai^tes tlie power ajvdt fuel flow of an internal cosabuetion piston 
engine. It can be used to determine the engine size required to meet an 
aircraft performance requirement/ or to predict the engine performance at a 
given operating point. It uses generalized dimensionless relationships between 
power and rim, auid between corrected power amd altitude so as to predict the 
full throttle power of an angina for any realistic combination of rptn and 
altitude. 

When operating losses are ignored, the power available from a piston 
engine is proportional to the product of displacement, rpm, and brake mean 
effective pressure or throttle setting. The losses increase with the rpm, 
so the power of a specific engine varies as shown in Figure IV. 1.4. The 
non-dimensional relationship between power and rpm container in PWRPLT is 
illustrate-*, in Figure 17.1.5. Also illustrated is the relationship between 
corrected power and altitude. 

Super charged (and turbo cheurged) engines maintain their rated sea level 
power up to some critical alt-^ ;ude above which maximum power decreases. 

Maximum power at altitudes above the critical altitude is related to maximum 
sea level power by 


where 


HPM 


(SIGMA - .117) 
(SIGCRT - .117) 


X HPMSLS 


SIGMA 


air density at altitude 
air density at sea level 


SIGCRT « 


air density at critical altitude 
air density at sea level 


as discussed in Principles of Aerodynamics, authored by Dwinnel (McGraw-Hill, 
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1949) 






FKSURE IV.1.5PART POWER AND ALTITUDE PERFORMANCE OF 
AIRCRAFT PISTON ENGINES 
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Til* op*rating point of a particular angina (fixad displacaaant) i* 

■pacified by the engine RPM and engine manifold pressure. In aircraft equipped 
with con8t^mt speed propellers, the pilot controls nmmifold pressure with 
the thrjttle; he controls engine speed by setting the propeller governor 
which adjusts the propeller blade angle such that the propeller absorbs the 
power developed by the engine at the desired throttle setting and engine RPM. 
Generally, either the engine operating point (power and RPM) or the aircraft 
flight condition (altitude and airspeed) is known and the other must be determined. 

The analysis of engines with fixed pitch propellers is somewhat more 
complicated. In this case the only power control is the tlirottle position. 

For a given throttle position and aircraft airspeed and altitude, the engine 
operates at that RPM at which power absorbed by the propeller equals the power 
produced at the crankshaft. The pilot may indirectly control RPM by adjusting 
the throttle position and aircraft airspeed (trim control) such that, at the 
desired engine RPM, the power developed by the engine is absorbed by the 
propeller. 

During engine sizing at cruise, the rated sea level horsepower of the 
engine is computed from the power required to drive the propeller (HPWR) and 
the input cruise engine power setting (PCPOWR * power at cruise/SLS power) : 

HPMSLS « HPWR/PCPOWR 

Reciprocating engine fuel consumption is expressed in PWRPLT as an empirical 
function of engine displacement and engine power setting. These relationships 
are illustrated for naturally aspirated and turbocharged anginas in Figures 
IV. 1.6 and IV 1.7, respectively. Fuel consumption at power settings less than 
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or ifequal to 78 par cant of aiaximum «aa leval powar is ccwputad ajssuming a 
lean fuel mixture. 

The input option KOBE (« 1 to 9) specifies whether the angina sise or 
the engine perfomumce is found, and the other arguments in the call state- 
ment are 

HPM “ maximum power available at altitude (output) 

HPMSLS ■ maximum sea leval power at full throttle, !'’?«ximimi rpm 
(output) 

HPWR » actual power at operating throttle and rpm (input or output) 

HPAVLB »» maximum full throttle power at altitude and at operating rpm (output) 
PCPOWR «■ per cent rated power (input or output) 

PCRPM » per cent maximum rptn (input) 

DELTA * operating static pressure ratio (input) 

RTHET - square root of operating static temperature ratio (input) 

H « altitude, ft (input) 

KSPCHG * supercharger flag (input) 

BSFC “ brake specific fuel consumption, lb per hr par hp (output) 

The first six of these quantities are illustrated in Pigtire IV. 1.4. 

IV. 1.2. 4 Turboprop Engine Performance, Subroutine TURBEG . Nearly half 
of this 270 card subroutine is numerical data, descriptive of the AIKESEARCH 
TPE331-1 turboprop engine. The program is used to scale this engine at a 
given flight condition and operating point, or to compute the performance of 
a given size engine at a specified operating point. The performance parameters 
of interest are shaft power, fuel flow and jet thrust. 
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Tha angina operating point and the flight condition togathar define the 
engine performance. Tha flight condition is specified by tha first three 
call parameters: 

PO » static praaaura, lb par sq ft (input) 

TO ■ static temperature, deg R (input) 

SUN » Mach number (input) 

Engine performance is measured by the next several call variables, i.a., 
HPWR, HPAVLB, HPM emd HPMSLS are as defined in the previous section while: 
PCNCR » per cent corrected maximum rotor speed (input or output) 

PCN » cent maximum rotor speed (output) 

T4SET « turbine inlet temperature, deg R (input) 

VTF « fuel flow, lb per hr (output) 

FN “jet thrust, lb (output) 

The remaining call parameters eire 

XNMAX “ maximum engine RPM (input) 

GR “ gear ratio; maximum propeller rpm/maxiraum engine rjaa (output) 
MODEP “ 0, cruise operation 

“ 1, takeoff operation (input) 

KODE “ 1 to 7, engine sizing options (input) 

The performance of the reference engine is scaled by straightforward 
means according to the value input for KODE: i.e., 

KODE “ 1 engine is being sized at a given flight condition, PCNCR is 
input and T4 is either input or a function of PCNCR 



KODE * 2 or 7 arkgina,sis« fixed, PCNCR ia input, T4 is input, or 
T4/T2 4-8 a function of PCNCR 

3 or 4 same ajii KOBE * 1, except required power is a fraction of 
the pqwer sized at input value of PCNCR 
5 or 6 angina size fixed, PCNCR is input, determine T4/T2 
so as to balance required and available power 
The corrected performance figures of the TPE 331-1 engine, at three flight 
Mach numbers (0., .25 and .50), are taibulated in TURBEG as functions of the 
engine operating point. The actual performance, for arbitrary altitude and 
Mach number, is found by interpolation or extrapolation of the tabulated data, 
and by applying the correction factors. 

A major portion of the subroutine is concerned with the scaling of 
several performance parameters from the reference engine performance data. 

The performance is specified by the horsepower, fuel flow, airflow and jet 
thrust, all of which vary linearly with engine size. Such operating variables 
as tvurbine inlet temperature and corrected rotor speed may be specified 
independently as constraints on the scaled engine. 

The performance of the TPE 331 engine is scaled by the ratio of the sea 
level static horsepower of the scaled engine (HPMSLS - determined during engine 
sizing or input) to the sea level static horsepower of the TPE 331 (HPSLRF) : 


horsepower; 


HPWR. 


SCALED 


»"«\PE 331 ^ 


HPMSLS . 

HPSLRF^ 


jet thrust; 

FN 

SCALED 

*« FN 

”tPE 331 

.HPMSLS 
^ ^HPSLRF 

fuel flow; 

WP 

SCALED 

m V7F 

TPE 331 

HPMSLS 

^ 'hpslrf 


) 

) 
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The acalad airflow ia used to computa tha maximtoa angina spaad of tha scaled 
angina frcai which tha scaled angina’s gaair ratio is determined i 


angina spaed: 


geax ratio; 


^%AXgcALED 

GR 

SCALED 


^%AXtPE331 


^\PS331 

WA 

SCAliED 


^Ve 331 Rpj^ 


TPE331 

•SCALED 


Note that maximum propeller shaft speed remains unchemged as the engine is 
scaled. 


IV. 1.3 Propeller Characteristics, Subroutine ENGDAT 
This subroutine controls the calling of the propeller related routines 
and deals with four aspects of the propeller requirements: 

1. Performance option - finds power/ thrust/blade angle relationship 
at given flight condition and propeller rpm 

2. Cost option - finds cost of propeller and gearbox 

1 3. Weight option - finds weight of propeller and gearbox 

4. Noise option - finds noise of propeller and gearbox 
The indicator KODE specifies which of the four options is desired, according 
to the value given to this argument; 

KODE « 1 to 10: performemce option 

11 to 20; cost option 
21 to 30: weight option 

31 to 40: noise option 

' ther input-output arguments of the program are, in order; 
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GR »* gear ratio, prop rpra/angina rpm (input) 

DROT * propeller diameter, ft (input) 

THRUST* thrust, lb (innut/output) 

SHP * shaft horsepower (input/output) 

EFPP » propeller efficiency (output) 

VKTS » airplane velocity, kts (input) 

RORO * ratio of air density to sea level density (input) 
lERROR* error indicator (output) 

ENP * number of engines (input) 

PO,TO » static pressure and temperature at altitude (input) 

AFX * activity factor (input) 

BLX * number of blades (input) 

The propeller performance option is the most ccmiplex. Whenever propeller 
performance is computed, the propeller geometry, RPM, zind aircraft airspeed 
are known. Together they define the advance ratio J; 

^ 101.4 X Airspeed (kts) 

" RPM X Diameter (ft) 

For a given advance ratio, the propeller performance problem can taka one of 
three forms: 

1. knowing thrust and advance ratio, find blade angle and power 

2. knowing power and advance ratio, find blade angle zind thrust 

3. knowing blade angle and advance ratio, find thrust and power 

The first problem is encountered by a constant speed (variable blade emgle) 
propeller during cruise at a specified altitude and Mach number. The second 
problem occurs with a constimt speed propeller during take-off, climb, and 
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cmi«e at a spacifiad power aatting. The firua.1 case is used to compute the 
performance of a fixed pitch propeller during all flight segments. In all three 
cases ENGDAT calls subroutine PERFM which consists of generalized propeller 
performance tables. 

IV. 1.3.1 Subroutine PERFM . PERFM is the propeller performance subroutine 
based on Haiailton Stamdard methods described in NASA CR-2066. For a given 
propeller geometry it relates power Jiind thrust coefficients f advance ratio, and 
blade emgle. Correction factors zire applied to account for differences in 
number of blades per propeller, activity factor per blade, and blade integrated 
design lift coefficient. 

Subroutine PERFM is nearly 500 cards in length, but about half of the program 
is numerical data, descriptive of propeller relationships, including blade 
geometry, propeller aerodynamics, power coefficients, etc. The remainder is 
concerned largely with the interpolation of this input data for the particular 
propeller input characteristics. Use is made of the utility subroutines 
BIQUAD and UNINT, for biquadratic interpolation of Y(X), and for uniform 
four-point interpolation. PERFM is called by the propeller ENGDAT. The input 
parameters to subroutine PERFM are 

IW * 1/ propeller power coefficient input 

« 2, propeller thrust coefficient input 
» 3, reverse thrust being calculated 
>■ 4, propeller blade emgle input 
ZJl « propeller advance ratio 

AFT *« activity factor per blade 

BLADT » number of blades 
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CLI ** prop«ll«r blada integratad dasign coafficiant 
ZMS(l) * forward Mach number 
ZMS(2) » tip Mach number 
CP if IW « 1, power coefficient 
CT if IW « 2, thrust coefficient 
BULLL if IW«4, blade angle 
Output parameters from subroutine PERFM are 
BLLLL and CT if IW » 1 
BLLLL and CP if IW - 2 
CT and CP if IW » 4 

ASTERK « error flag if there is problem calculating propeller 

performamce 

XFT * ccxnpressibility correction factor 

There are several numerical tests in the code related to being within the 
limits of the tabular data; i.e., if the error flag LIMIT is returned as 1 
by either BIQUAD and UNINT, the data is outside the lower end of the tabular 
data. If LIMIT returned as 2, the data is outside the high end of the tabuleir 
data. An error message indicates in what table the problem occurs. 

IV. 1.3. 2 Propeller Costs, Subroutine COST . This subroutine estimates 
propeller costs according to 1970 or 1980 manufacturing technologies and is 
called by the propeller ENGDAT , The following parameters are the input; 

WTCON » type of propeller; « 1, fixed pitch propeller 

» 2, constant speed propeller 
« 3, constant speed, full feathering propeller 
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4, conatzmt sps®d|. full fsatharixig/ 


daioing propallar 

* 5, consteuit speed, full feathering, 
deicing, reversing propeller 

BLADT *« number of blades 
CLFl » same as XCLFl 

CLP “ same as XCLP 

CK70 « same as XCK70 

CK80 • same as XCK80 

CAMT » nvunber of propellers produced (optional) 

WT70 « propeller weight, 1970 technology, lb 

WT80 « propeller weight, 1980 technology, lb 

lENT « 1, initialization for propeller cost factors 
» 2, computes propeller cost 

The output quantities are found as numerical functions of these* 

CQUANCD* number of propellers produced, 1970 technology default is a 
function of propeller type or can be input as CAMT 
CQUAN(2)» number of propellers produced, 1980 technclogy default is a 
function of propeller type or can be input as CAMT 
COST70 » propeller cost, 1970 technology 
C0ST80 « propeller cost, 1930 technology 

CCLFl “ learning curve factor ~ default function of propeller type or 
can be input as XCLFl 

CCLF «* leeurning curve factor - default function of propeller type or 
can be input as XCLF 
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?omuXa« tut«d ar« 


0 75 

COST70 or COST80 »« C « ZJT (3B +2) 

D.t.uU V«lu. of CCK70 or CCK80 - C, . f ♦ B) 

whare 

COST70 or COST80 « C « Average O.E.M. propeller cost for a number of 
units/year r $/lb 

CCK70 or CCK80 “ ■ single vinlt O.E.M. propeller cost/ $/lb. 

Z « ^/^l 

CCLP « LF » learning curve factor for a number of units/year (default 

« 1 . 02 ) 

CCLPl » learning curve factor for a single unit (default » 3.2178) 
BLADT « B ■ number of blades 

F « single unit cost factor 

E » empirical factor 

(NOTE: reference Figure IV. 1.8 for LF zmd LFj_ values baaed on an 89 per cent 

slope learning curve.) 

Constants used in these equations are 


WTCON or 
NTYP 

1970 

CQUANd, VJTCON) 

1980 

CQUAN(2, WTCON) 

P 

E 

Quantity 

F 

E 

Quantity 

1 

3.5 

1.0 

1910 

3.5 

1.0 

2230 

2 

3.7 

1.5 

2810 

3.7 

1.5 

5470 

3 

3.2 

3.5 

1030 

3.2 

3.5 

1990 

4 

2.6 

3.5 

295 

3.5 

3.5 

680 

5 

2.0 

3.5 

65 

3.4 

3.5 

368 
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IV. 1.3. 3 G earbox Weight, Coat anid Noi«« Ch^lract«riatic«, Subroutine 
GEiVRBX . The gearbox weight, cost eind noise can all be found by this subroutine 
based on the 1973 Hamilton-Standard study, NASA CR-114665, "Q-FAN for General 
Aviation." The last of the input arguments is the flag MODE, which takes the 
values 1, 2, or 3, respectively, according to whether weight, cost or noise 
is to be computed. Other input parameters are 
XNMAX - maximum engine rpnv 
PCRPM ■ fraction of maximum rpm 
SHP » shaft horsepower 
DHOT ■ diameter of propeller, ft 
i'.GR - gear ratio (propeller rpm/engine rpm) 

CATN « aircraft type 
KWKITX- write flag 

The gearbox parameters then follow as numerical functions of these quantities 
iu»d are output as 

GPNDB * maximum gearbox noise at 500 ft, decibels, PNdb 
GDBA « maximum gearbox noise at 500 ft, decibels, DBA 
WTGB * gearbox weight (including mount and afterbody) , lb 
CSTGB « gearbox cost, $ 

It may be noted that the noise is proportional to a quantity X " 10 log (SHP) 
which is measured in units of decibels. The scale factor is 10 instead of 20 
because the power is p>rojK)rtional to the square of the noise, which introduces 
a factor of 2 into the decibel representation of the noise. 
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The gearbox ascembly incltidss 

• Housing 

• Bearings 

• Planetary Gearing 

• Tailshaft 

• Afterbody 

• Lube and scavenge punip (single or 2-stage gearing as required) 

• Fan accessory Drives 
Weights are given by 


Single-stage; 


W_ 


8 


nO.84 


(SHP)P 

TS 


+ 0.6 

AFTERBODY 


(for gear ratio > .20) 


Two-Stage: 


W. 


T 


10.6 


(SHP)D 

TS 


0.84 


+ 0.6 D^ 

AFTERBODY 


(for gear ratio .20) 


Gearbox coat is given by 


CSTGB » * Z * WTGA + 13.5 * Z * WAFTB 

where 

« first unit cost ($/lb) 

“ 150 for single-stage planetary 
■« 180 for two-stage planetary 
and 

Z - (LF/LFj^) 
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wh«r« LP « learning curve factor for nvtmbar unit«/yr 
» learning curve factor for first unit 
Typical values for Z are 


239 

NTYP - 

1 and 2 

283 

m 

3 

338 

M 

4 

374 

m 

5 


Noise levels are predicted by 

Single stage planetary: GPNDB ■» 31.0 + 10 log SHP 
Two-stage planetary: GPNDB * 34.0 + 10 log SHP 

and • GDBA « GPNDB - 11 

IV. 1.3. 4 Propeller Weights, Subroutine WAIT . Propeller weight is 
estiniated in this 30-card subroutine as a numerical function of seven input 
parameters : 

NTYP ■ IWTCON « airplane propeller type (1 to 5) 

ZMWT « Mach number correction to propeller weight 

BHP « breike horsepower 

DIA ■ propeller diameter/ ft 

AFT ■ activity factor per blade 

BLADT- number of blades 

TIPSPD- tip speed, ft per sec 

Then, according to several straightfoxrward but nonlinear functions, the 
output parameters are simply: 

vrrVO *• propeller weight, 1970 technology, lb 
WTRO *• propeller weight, 1980 technology, lb 
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Equations amployad ara 



where 

WT70 or WT80 « « propeller wet weight, lbs. (excludes spinner, de- 

icing euid governor) 

DIA » D « propeller diameter, ft 
BLADT « B « number of blades 
AFT * A.F. »• blade activity factor 

VmTp 

N » propeller speed, rpm (take-off » — j * TIPSPD 

BHP » SHP » shaft horsepower, HP (take-off) 

ZMWT « M » Mach number (design condition; maximum power cruise) 
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Propallar types asscxjiated witti above and are as follows: 

(1) all fixed-pitch props 

(2) McCauley non-counterweighted, non-feathering, constant spaed props 

(3) All Hartzell, all Hamilton Staindard small props, and feathering McCauley 

(4) Fiberglass-bladed, constant speed, countarweighted, full feathered 

(5) Fiberglass-bladed, constant-speed, double-acting (non-counterweighted), 
full feathered, reverse 

IV.1.3.5 Propeller Noises, Subroutine ZNOISB . Most of this subroutine is 
numerical data defining the noise generated by the propeller. The subroutine 
has eight other input quantities, and they are, in order; 

BLADT « number of blades 

DIA ■ propeller diameter, ft 

TIPSPD * propeller tip speed, ft per sec 

VKTAS » aircraft velocity, kts 

BHP « brake horsepower of engine 

DIST “ slant distance to observer, ft 

* "^"^sea level'^^ambient 
XNOE » number of engines 

These quantities are then used to develop the output, which is 
SPL » total perceived noise level, PNdB 
SPLX » total sound pressure level, DBA 
Equations employed are 

SPL - XLl + XL2 + XL3 + XL4 + PNLD 

where 
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XLl »« ref level from Figure IV. 1,9 at raf condition 

XL2 » diameter and blade correction » «20[log DIA/10,5 + log BLaD/4,0] 

XL3 »« distance ev .rrection « -20 log DIST/500. 

XL4 « number of engines correction « 10 log XNOE 

PNLD ■ perceived noise adjustment - table look-up function of BLADT, 

DIA, and helical tip Mach number 

IV. 1.3. 6 Propeller Driven Aircraft Noise Control Routine, Subroutine 
PNOYS . PNOYS controls noise computation for propeller aircraft. This routine 
is called by MAIN to get noise. As shown in Table IV. 1.1 this is the third 
principal subroutine to gat propeller 2 ind gearbox noise,; engine noise is 
computed with a call to ZNENG which will be described in the next subsection. 
PNOYS is about 45 statranents in length, eind it is directed principally by the 
argument KNOYS: i.e., 

KNOYS « 0, both Mach number and altitude are given, and noise is determined 

KNOYS *« 1, altitude and power setting given, and Mach numiber and 
noise are determined. 

In the latter case, Mach number is found by a call to subroutine ASPEED to 
find Mach number capability at altitude and power setting, and the computations 
after this point are independent of KNOYS. 

Most of the mechanical and operational input parameters are familiar 
from other subroutines; e.g., DPROP is the propeller diaimeter, and GRATIO is 
the gear ratio, etc. Other output from PNOYS is returned by subroutines ENGINE, 
GEAJRBX and XNENG, the last of which deals with noise characteristics, and it is 
described in the next subsection. 
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FIGW^ IV. 1.9 - BASIC NOISE LEVEL CURVE FOR 4"BLADED^ 10.5 FT DIA. PROPELLER 9 500 FT. 















IV. 1.3. 7 Bnglna Noise Charactarlgtlcg ^ Subroutine ZNBNG , This subroutine 
computes the noise characteristics of the following types of engines 
JTYPE « 1, piston engine 

* 2, water cooled rotary ccajibustion engine 
» 3, turboshaft engine 

The noise at 500 feet distance is found assuming that the aircraft Mach number 
is .1, and on the basis of experimentally derived numerical expressions illus~ 
trated in Pigxires IV, 1.10 to IV. 1.12. Other descriptive inputs to the subroutine 
are 

SHP * engine shaft horsepower 

XNMAX « maximum engine rpm 

PCRPM » ratio of operating rpm to maximum rpm 

NOE * nimiber of engines 

GGR * gear ratio (prop rpm/engine rpm) 

The output of the subroutine are the noise levels 
EPNDB »* perceived noise level, PNdB 

EDBA » weighted level as measured on the A-scale, for which the noise 
levels are reduced at low and high frequencies, dbA 
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FIGL-RE IV. 1. 10 - NOISE OF UNMUFFLED GAS TURBINE ENGINE 




30 60 100 300 600 1000 

KW 

'engine shaft power 


PTRIIRP TV T TT - NOISE OF UNMUFFLED WATER COOLED ROTARY 
i-lbUKt lY.l.il COMBUSTION ENGINES 
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FIRING FREQUENCY 


#N^*NCYL 


WHERE: N ^ 0,5 FOR FOUR CYLINDER 
N “ 1,0 FOR TM) CYLINDER 

NCYL ® NUMBER OF CYLINDERS 


RPM 

NCYL 

F.F. 

m} 

6 

170 

2700 

6 

135 

1500 

6 

75 

2700 

4 

90 

2010 

4 

6G.7 


dBA - PNdB - 12. 


FIRING FREQUENCY 


MAX! 

500 FT 
SIDELINE 
PERCE 
NOISE 
LEVEL 
IN PNdB 


EXHAUST 

NOISE 




INLET 

NOISE 


60 100 SHP 300 

60 100 300 600 

kw 

ENGINE POWER 


600 1000 ^2000 
1000 2000 


FIGURE IV. 1. 10 “ NOISE OF UNMUFFLED PISTON ENGINES 
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GASP - GENERAL AVlATfON SYNTHESIS PROGRAM 


VOLUME IV - PJE^OPULSiCON 
PART 2 - USER'S MANUAL 


JANUARY 1978 


Pr«par«d for 

miCJKAL ASEONAUTICS AND SPACE ADHTNISTRATrCSf 
Ames Research Center 
Moffett Field, Califomi* 

Under 

CONTRACT NAS 2-9352 


AEROPHYSICS RESEARCH CORPORATION 



IV. 2 PROPULSIOSr MODEL USER'S MANUAL 


Th« propulaion niod#! aubroutinaa are very maaeroua/ whether the syateia 
ia of turbofan or propeller in form. The present section alphabetically 
tabulates and defines the input/output parameters for all the propulsion 
subroutines, with turbofan programs followed by propeller programs. These 
tabulations follow the order given below. 

Turbofan Svibroutinea Propeller Subroutines 


ENGDTl-7 

COST 

PNOYS 

ENGDTT 

ENGDAT 

PWRPLT 

ENGINE 

ENGINE 

TURBEG 

ENGSZ 

ENGSZ 

WAIT 

NACDG 

GEARBX 

ZNENG 


PERFM 

ZNOISE 


Many of the subroutines are devoted principally to the tabulation of a 
specific propulsion system performance data, and in these cases only a few 
additional input-output queintities eire listed. On the other hand, the very 
long programs (ENGSZ and PERFM, for example) eire associated with many such 
parameters. As will be seen, the propeller systems generally require a greater 
number of input peurameters them do the turbofem systems. 

The seven subroutines listed below are not tabulated because they deal 
with numerical perforraemce cheiracteristics of a group of aircraft turbofan/ 
turbojet jet engines; 
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ENGDTl * GE CJ610-6 Turbojat 
ENGDT2 * Garrett TFE 731-2 Turbofan 
ENGDT3 » UACL JT15I>-1 Turbofan 
ENGDT4 » Lycoming ALF-502 Turbofan 
ENGDT5 « GE CP34 Turbofan 
ENGDT6 » GE TF34 Turbofan 
ENGDT7 • QCGAT Turbofan 
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FIGURE IV. 2.1 SUBROUTINE ENGDTT (TURBOFAN) —INPUT 


VARIABLE 

DESCRIPTION 

FTUKOT 

jKJwer sottiiv) as a fraction of maximum 

HN 

altitude/ ft 

KENG 

engine power setting indicator (0 to 7) 

NALT 

number of altitudes in tables 

NMN 

number of Mach numbers in table 

NT4 

number of turbine inlet temperature ratios T4/T2 in table 

P2 

total pressure, lb per sq ft 

SFNIDL 

idle specific thrust, lb per Ib/sac 

1 SM 

Mach number 

T2 

total temperature, deg R 

T4MAX 

maximum turbine inlet tmperatura , dag R 

T4MC 

turbine inlet temperature in cruise configuration, deg R 

T4MCL 

turbine inlet temperature in climb configuration, deg R 

WAMAP 

maximum sea level static airflow of reference engine at 


100 per cent corrected rotor speed, lb per sec. 
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FIGURE IV. 2. 2 SUBROUTINE ENGDTT--OUTPUT 


VARIABLE 


DESCRIPTION 


fuel air ratio 


PCWAC 


thrust, lb 

ratio of corrected airflow to maximum sea level static 


airflow 

specific fuel consumption, lb per hr per lb of thrust 
specific thrust; lb per lb per sac of airflow 
airflow, lb per sec 


fuel flow, lb per hr 
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7IGUKE IV. 2, 3 SUBROUTINE ENGINE (TURBOEAN) —INPUT 


VARIABLE 

DESCRIPTION 

FAR 

fuel-air ratio 

FN 

thrust, lb 

IIN 

altitude, ft 

lENGSC 

engine cycle indicator 

ISEGX 

mission segment indicator 

KENG 

engine power setting indicator 

KNAC 

nacelle drag/sizing indicator 

KODEAC 

acceleration segment power setting indicator 

KODECL 

climb segment power setting indicator 

KODETO 

takeoff segment power setting indicator 

KODETR 

turn segment power setting indicator 

KWRITE 

write indicator 

PCWAC 

ratio of corrected airflow to meocimum sea level static 


airflow 

PR 

inlet pressure recovery factor 

PO 

static pressure, lb per sq ft 

SFC 

specific fuel consumption, lb per hr per lb of thrust 

SMN 

Mach nvunber 

TO 

static temperature, deg R 

WASLS 

sea level static airflow, lb per sec 

WG 

design gross weight, lb 


wing loading, lb per sq ft 
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FIGURE IV. 2. 4 SUBRjOUTINB ENGINE— OUTPUT 










FIGURE IV. 2. 5 SUBROUTINE ENGSZ (TURBOFANji— -INPUT 


VARIABLE 

DESCRIPTION 

CDNI 

nacelle skin friction coefficient 

CLMXLD 

maximum lift coefficient in landing configuration 

CLMXTO 

maucimum lift coefficient in takeoff configuration 

DCLTO 

lift coefficient increment in takeoff configuration 

DRG 

total drag , lb 

DRGNCl 

nacelle drag, lb 

EM 

Mach number 

EMTURN 

Mach number in steady turn 

ENP 

number of engines 

H 

altitude, ft 

HBTP 

hub to tip ratio or fan 

HPORT 

airport altitude, ft 

HTURN 

turn altitude, ft 

ICRU 

indicator used when KNAG « 0 (nacelle drag « thrust loss) 

lEGWGT 

indicator to determine if engine dimensions are to be 


calculated by ENGWGT 

JENGSZ 

engine sizing option indicator 

KNAG 

nacelle drag sizing indicator 

NACDRG 

indicates if engine dimensions have been calculated 


external to ENGSZ 

NPC 

computation indicator 

NTYE 

type of engine indicator 

NTYP 

type of propeller indicator 
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VARIABLE 

DESCRIPTION 

PCWAC 

ratio of corrected airflow to maximum SLS airflow 

PO 

static pressure, lb per sq ft 

RCCRU 

required cruise rate of climb capability, ft per min 

RWCRTO 

ratio of cruise weight to takeoff weight 

SFNSLS 

specific thrust, sea level static, lb per lb par sec 

SMID 

fan face Mach number, assuming one-dimensional flow 

SW 

wing area, sq ft 

TDELTO 

temperature deviation from standard day, for takeoff 


engine sizing 

THIN 

input thrust per engine (JENGSZ « 4) 

TO 

static temperature, deg R 

WF 

fuel flow, lb per hr I 

WG 

gross weight, lb 

WGS 

wing loading, lb per sq ft 

WTRFAC 

weight during turn divided by maximum gross weight 

XCKN 

nacelle form factor 

XLFTRN 

turn load factor 

XL^DE 

nacelle length to diameter ratio 

XTO 

actual takeoff distance, ft 

XrORQ 

required takeoff distance, ft 
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FIGURE IV. 2. 6 SUBROUTINE ENGSZ-— OUTPUT 


VARIABLE 

DESCRIPTION 

AE 

flow area at fan face, sq ft 

ANAC 

nacelle surface area, sq ft 

JDNAC 

nacelle drag coefficient 

>JLTLMT 

limit lift coefficient during turn 


Mach ntirnber 

i;np 

number of engines 

KNSLS 

sea level static thrust, lb 

rue 

control flag 

IPART 

FAR regulation part indicator for climb sizing 

ISEGX 

flight segment indicator 

'i 

-ITRS25 

1 engine sizing for turn indicator 

r.SIZE 

1 

engine sizing flag for takeoff 

liELl 

Reynolds number per foot 


current aircraft weight, lb 

.VASLS 

sea level static airflow, lb per sec 

XLN 

nacelle length, ft 

YCLB 

used in subroutine NACDG for including nacelle drag as 


engine thrust loss 

YDRG 

used in subroutine NACDG for including nacelle drag as 

: 

engine thrust loss 

ZLQD 

lift to drag ratio 
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FIGUBE IV. 2. 7 SUBROOTIMS HmCEKS (TTOBOrAN) — INPOT 


VARIABLE 

DESCRIPTION 

AN AC 

surface area of nacelle, sq ft 

CDNl 

reference drag coefficient of nacelle 

ENPP 

number of engines 

ICRU 

engine siaing/performanca flag 

KWRITE 

print indicator 

PO 

static pressure, lb par sq ft 

RELI 

Reynolds number per unit length, per ft 

SMN 

Mach number 

SWING 

wing uiea, sq ft 

XCKN 

nacelle form factor 

XLN 

nacelle length 

XLQDE 

engine (nacelle) length/diamater ratio 

YCLB 

thrust increment required to provide cruise climb margin, 
lb. 

YDRG 

cruise drag, level flight, lb. 
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ffIGUKB IV. 2. S SUBBOOTINS NACDG — CXTOPOT 


VARIABLE 

DRGNCi 

SFC 

SFN 

WA 


DESCRIPTION 

drag of one nacelle/ lb. 

specific fuel consumption, corrected for nacelle drag 
specific thrust Ib/lb/sac FUA^a., corrected for nacelle 
drag 

airflow, Ib/sec, corrected for nacelle drag 






FIC3URE IV, 2. 9 SUBa3OTIKE EMGSZ (FMKmiJSE) -~IWOT 








VARIABLE 


DESCRIPTION 


NTYP 

propeller type indicator 

PCPCR 

per cent maximuro sxjwer in cruise for reciprocating 


engines 

PCRCR 

per cent maximum rpoi in cruise for reciprocating engines 

PR 

nacelle inlet pressure recovery factor 

PO 

static pressure, lb per sq ft 

RCCRU 

required cruise rate of climb capability, ft par min 

RELI 

Reynolds number per imit length, per ft 

RPM 

engine speed, rev per min 

RWCRTO 

ratio of weiaht at cruise to gross weight 

SW 

wing area, sq ft 

TDELTO 

temperature deviation from standard, deg F 

TSPDMX 

maximima propeller static tip speed, ft per sec 

TO 

Static temperature, deg R 

WG 

gross weight, lb 

WGS 

wing loading, lb per sq ft 

XCKN 

nacelle form factor 

XLQDE 

nacelle length to diameter ratio 

XNMAX 

rticiximum engine speed, rpm 

XTO 

actual takeoff distance, ft 
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PIGUIE IV. 2. 10 SUBWOTIKB WGSZ'-KXJTFOT 


VARIABLE 

' " 

DESCRIPTION 

ANAC 

nacelle surface area, sq ft 

CDNAC 

drag coefficient of nacelles based on wing area 

DPROP 

propeller diameter, ft 

EM 

Mach number 

GRATIO 

gear ratio, propeller rpn/engiixa rpm 

H 

altitude, ft 

HNCRU 

cruise altitude, ft 

HPMSLS 

sea level static maximum horsepower 

IDC 

special purpose indicator 

I PART 

FAR part 23/25 climb requirement indicator 

ISEGX 

segment iMicator 

KFPTCH 

fixed pitch propeller indicator 

ICSIZE 

engine takeoff sizing indicator 

PCPOWR 

fraction of maximum power 

PCRPM 

fraction of maximum rpm 

PO 

static pressure, lb per aq ft 

RELI 

Reynolds number per unit length, per ft 

TO 

static temperature, deg R 

XliN 

nacelle length, ft 

XTORQ 

takeoff distamce required to clear 35 ft altitude, ft 

XLQD 

lift to drag ratio 
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WX&jm IV. 2.11 SUBROUTINE ENGINE (PROPEIIJ®) —INPUT 


VARIABLE 

DESCRIPTION 

AF 

propeller blade activity factor per blade 

BL 

number of propeller blades 

BLANG 

propeller blade angle at 3/4 radius, deg. 

COD 

QFAN shroud fineness ratio 

CP 

propeller power coefficient 

CPRDP 

propeller cost, $ 

Q2FT 

cost of "Q-FAN" propulsor, $ 

CTI 

initial estimate of propeller thrust coefficient 

DPRDP 

propeller diameter, ft. 

ENP 

number of engines 

FT 

propeller slipsti '• or» fraction of thrust lost 

GR 

(T /T ) 

effective isolated 

gec.r ratio; propeller rptn/engine rpm 

H 

altitude, ft. » 

HCRIT 

critical altitude for turbocharged engine 

HPAVLB 

horsepower available 

HPM 

maximum horsepower at given altitxide 

HPMSLS 

meiximum horsepower at sea level standard conditions 

ISEGX 

mission segment indicator 

KENG 

power setting indicator 

KODE 

engine sizing, power and flight condition options, defined 


in ENGINE 

KODECR 

reciprocating/turboprop engine cruise sizing option; see 
namelist INPROP 
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VARIABLE 

DESCRIPTION 

KUDKTH 

rociprocatinq/turbopi'op enqino throttling options 

KSPCHG 

supercharger indicator; see namelist INPROP 

KWRITE 

subroutine print option; see namelist INGASP 

NTYE 

engine typ>e indicator; see namelist INGASP 

NTYP 

propeller typ>e indicator; see namelist INPROP 

PCNCCL 

per cent corrected rotor speed at climb, turboprop; see 


namelist INPROP 

PCNCCR 

per cent corrected rotor speed at cruise, turboprop; see 


namelist INPROP 

PCNCTO 

per cent corrected rotor speed at takeoff, turboprop; see 


namelist INPROP 

PCPCL 

per cent maxinuan power in climb, rociproca^ting engine; sec 


luunelist INPROP 

PCPTO 

per cent maximum power at takeoff, reciprocating engine; 


see namelist INPROP I 

PCRCL 

i-.-r cent miiximum rpm in climb, reciprocating engine; see 

M 


namelist INPROP 

PCRCR 

per cent maximum rpm in cruise, reciprocating engine; see 


namelist INPROP 

PCRTO 

per cent maximum rpm at takeoff, reciprocating engine, sei* 


namelist INPROP 

PO 

static pressure, lb par sq ft 

SHP 

shaft horsepower required to turn propeller 
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VARIABLE 

SMN 

THRUST 

TO 

TSPDMX 

T4STCL 

T4STCR 

T4STTO 

WQFT 

XJ 

XNHAX 


DESCRIPTION 

Mach number 
thrust# lb 

standard atmospheric temperature, deg R 
maximum propeller tipspeed, ft per sec 
turboprop turbine inlet temperature at climb, deg R 
turboprop turbine inlet temperature at cruise/ deg R 
turboprop turbine inlet temperature at takeoff, deg R 
weight of "Q-FAN" propulaor 
advance ratio 
i maximum engine spaed, rpm 


IV-2 

16 




FIGURE IV. 2. 12 SUBROOTINE ENGINE — OUTOOT 


VARIABLE 


DESCRIPTION 


propeller thrust coefficient 


HPWR 


engine power output/ hp 


KERROR 


error indicator 


KFPTCH 


fixed pitch propeller indicator set by NTYP 


KWRITY 


write indicator 


PCPOWR 


engine power output, hp 


PCRPM 


percent of maximum engine rpm 


engine rpm 


SIGCRT 


density ratio at critical altitude for supercharged 


reciprocating engine 


THRUST 


l>ropeller thrust, lb 


TSFC 


thrust specific fuel consumption, lb par hr per lb 


fuel flow, lb per hr 
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FIGUKE XV. 2. 13 SUBKOUTINE BNGDAT (PIOPSLLKR) -—INPOT 


VARIABLE 

DESCRIPTION 

AFX 

propeller blade activity factor 

ASTERK 

error flag if there is problem in calculating propeller 
performance 

BEX 

number of propeller blades 

CAMT 

production quantity of propellers to be used (default or 
input in namelist INPROP) 

CLI 

design integrated lift coefficient of propeller 

DIST 

sl^mt diBt^mce to observer for noise 

DROT 

propeller diameter 

EM 

Mach number 

EMCRU 

! 

caruise Mach number 

ENP 

number of engines 

GR 

gear ratio, ratio of propeller rpra to angina rpta 

I DATE 

propeller technology level (1970 or 1980) 

KNAC 

nacelle drag indicator 

KODE 

engine performance options defined in subroutine ENGINE 

KWRITE 

print indicator 

PCRPM 

fraction of maiximum rpm 

RORO 

ratio of air density to standard sea level density 

TO 

temperat'jure, deg R 

VKTS 

airspeed, kts 

WKPFAC 

propeller weight adjustment factor 

XCK70 

single unit propeller cost (1970),$ per lb | 
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VARIABLE 


DESCRIPTION 


XCK80 

XCLF 

XCLFl 

XFT 

XNHAX 


single vinit propeller cost (1980) , $ per lb 

learning curve factor for yearly units (1.02) 

learning curve factor for single unit (3.2178) 

propeller compressibility correction (0, no compressibility 

maximum engine speed, rev per min 




FIGUKB rV.2.14 SOTKJOTINE 


VARIABLE 

DESCRIPTION 

AF 

i>ropeller activity 

BL 

mmtber of propeller blades 

BLAKG 

propeller blade emgle at 3/4 radius, dag 

BLLL 

propeller blade angle at 3/4 radius, dag 

CP 

power coefficient 

CPROP 

propeller cost,_ $ 

CSTGB 

gear box cost, $ 

CT 

thrust coefficient 

EFFP 

propeller efficiency 

I ERROR 

error indicator 

NTYP 

propeller type indicator 

SHP 

shaft horsepower 

THRUST 

total propeller thrust, lb 

WROPOl 

weight of one propeller, lb 

WTGB 

geau: box weight, lb 

ZJI 

advance ratio 

ZMWT 

cruise Mach number 
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FIGURE IV. 2. 16 SUBROUTINE O>ST-“0UTFUT 


VARIABLE 


DESCRIPTION 


CCK70 


average original equipment manufacturer 1970 propeller 


cost/ S per lb 


CCK80 


average original equipment manufacturer 1980 propeller 


cost/ $ per lb 


CCLF 


learning curve factor for a number of unite per year 


cx:lfi 


learning curve factor for a single unit 


COST70 


propeller cost, 1970 technology, $ 


COST80 


propeller cost, 1980 technology, $ 


CQUAN(I) 


number of propellers produced in year I 
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PICft3KE IV. 2. 17 SUBIOUTINE GEARBX (PIOPBUSK) -iNPOT 
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FIGURE IV. 2. 18 SUBROUTINE GEftRBX— CWTPUT 


VARIABLE DESCRIPTION 

CSTGB cost of gearbox, dollars 

GDBA gear box noise, dBA 

GPNDB gearbox noise , PNdB 

WTGB gearbox weight, lb 







FIGURE IV. 2. 19 SUBROUTINE PE.RFH (PROPELLER) —INPUT 


VARIABLE 


DESCRIPTION 


BLADT 


LIMIT 


activity factor per blade 
nvunber of blades per propeller 
blade angle, deg 

propeller blade integrated design coefficient 

thrust coefficient 

type of computation flag 

error return flag 

propeller compressibility factor 

advamce ratio 

propeller Mach number 
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FIGURE IV. 2. 20 SUBROUTINE PEEFM — OUTPUT 


VARIABLE 

DESCRIPTION 

AFCPE 

activity factor adjustment on effective power coefficient 

AFCTE 

activity factor adjustment on effective thrust coefficient 

ASTERK 

error retvu:n flag 

BLLL 

blade angle, deg 

CP 

power coefficient 

CPE 

effective power coefficient 

CTE 

effective thrust coefficient 

XFT 

propeller compressibility factor 
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FIGURE IV. 2. 21 SUBROUTINE PNOYS (PROPEUUER) —INPUT 


VARIABLE 

' 

DESCRIPTION 

□PROP 

propeller diameter, ft 

EM 

Mach number 

ENP 

number of engines 

GRATIO 

gear ratio; propeller rpm to angina rpm 

H 

altitude, ft 

HPMSLS 

maximum sea level static horsepower 

KNOYS 

noise calculation flag 

KWRITE 

output print flag 

NTYE 

engine type indicator 

NTYP 

propeller type indicator 

WF 

fuel flow, lb per hr 

WG 

aircraft gross weight, lb I 

XNMAX 

maximum engine speed, rpm I 







FIGURE IV. 2. 22 SUBROUTINE PNOYS—OUTPOT 








FIGURE IV. 2. 23 SUBROUTINE PWRPLT (PROPELLER) —INPUT 
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FIGURE IV. 2. 24 SUBROUTINE FMRPLT — OUTPUT 


VARIABLE 

BSFC 

HPAVLB 

HPM 

HPMSLS 

HPWR 

PCPOWR 

SIGCRT 


DESCRIPTION 

specific fuel consumption, lb per hr per hp 

loaxljnujn full throttle horsepower available at altitude at 
operating rpm 

ma"iinum horsepower available at altitude 
maximum sea level static horsepower 
engine power output, hp 
per cent maximum power 

critical air density ratio for supercharged engine 
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FIGURE IV. 2. 26 SUBROUTINE TURBEG — OUTPUT 










FIGORE ^.2.27 SUBROUTIHE WAIT (PROPELLER) —INPUT 


VARIABLE 

DESCRIPTION 

AFT 

activity factor per blade 

BHP 

braike horsepower 

BLADT 

number of propeller blajdas 

DIA 

propeller diameter, ft 

TIPSPD 

propeller tip speed, ft par sac 

WTCON 

parameter defining aircraft category 

2MWT 

Mach number correction on propeller weight 
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riGURB IV. 2. 28 SUBRCOTINB WAIT—OOTPUT 


VARIABLE 

WT70 

VJT80 



propeller weight, 1970 technology, lb 
propeller weight, 1980 technology, lb 








FIGURE IV. 2, 29 SUBROUTINE ZNENG (PROPELLER) —INPUT 


VARIABLE 



ESC 



ITYPE 


type of engine indicator set according to ntye 


number of engines 


PCRPM 


fraction of maximum rp© 


»aft horsepower 


XNMAX 


maximum engine speed, rev par min 
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FIGURE IV. 2. 30 SUBRO UT INE ZNSHG — OUTPUT 

VARIABLE DESCRIPTION 

EDBA noise level, dBA 

EPNDB I noise level, PNdB 







FIGURE IV. 2. 31 SUBROUTINE ZNOISE (PRDPELIER) —INROT 


VARIABLE 


DESCRIPTION 


brake horsepower 


BLADT 


number of propeller blades 


propeller diameter, ft 


DIST 


slant distance to observer, ft 


square root of taaparature ratio 


TIPSPD 


propeller tip speed, ft per sac 


VKTAS 


TAS, kts 


number of engines 
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FIGURE IV. 2. 32 SUBROOTIHE 22IOISE— OUTPUT 
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IV. 3 PROPULSION MODSL AND PROGRAMMER'S MANUAL 


Th® flow chart* of tha propulsion aubroutinas ara shown in alphabatic 
ordaTi. for both turbofan and propeller systatos. A total of 23 subroutines 
have been drawn and are presented in tha following alphabetic order: 


Turbofan Subroutines 


ENGDTl-7 

ENGDTT 

• ENGINE 

• ENGSZ 

• NACDG 


TYPIFIED BY ENGDTI 


Propeller Subroutines 


# 

COST 

• 

ENGDAT 

• 

ENGINE 

• 

ENGSZ 

# 

GEARBX 

• 

PERFM 

# 

PNOYS 


PWRPLT 

• 

TURBEG 

« 

WAIT 

« 

ZNENG 


ZN0I3E 


Tho flow charts do not include tha prasantation of numerical performance 
data, which is a l^lrga portion of many of the routines. On the other hand, 
the interdependence among the routines can be appreciated by nctiny the 
subroutines of each. 


lV-3 
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rv.3.1 Tiorbojat and Turbofan Routinas 

IV. 3. 1.1 Subroutlna ENGINE, Turbo j at aiva Turbofan Bnglna Parforaance 
at Specified Flight Condition . This routine cociputas turbojet and turhofan 
Lhrust, airflow, fuel flow, specific thrust, per cent corrected airflow, 
and thrust specific fuel consumption at a specified altitude aund Jiach number. 
Figure IV. 3.1 provides a detailed flow chart of this subroutine. Routine 
operation is described in Section IV. 1.1.1. 

The various engine sizing options are controlled by the indicator KENG. 
Engine type selection is controlled by the indicator lENGSC which taOtes one 
of eight values to call the engine type desired. Basic engine types avail- 
able are contained in the subroutines ENGDTl to ENGDT7 and ENGDTT. Subroutine 
NACDG is used to compute nacelle drag losses Engine iterations are carried 
out with the aid of subroutine ITRMHW, described in Section I. 1.3. 8. 
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ENGINE 



FIGURE IV. 3.1 - DETAILED FLOVTCHART, SUBROUTINE ENGINE 
^ 

3 











ENGIHE 


2 

















ENGINE 


S 
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rv.3.1,2 Subroutine JSWGSZ, Turbojat and Turbofan Sraglna 3ising . This 
routine cxmtrols turbojet end turbofan engine sising cotaputationa . Routine 
ftmctlon is described in Section IV. I. 1. Engine characteristics are deter- 
mined by calls to subroutine ENGINE. Nacelle characteristics are controlled 
by the indicator KNAC. Siting method is controlled by the indicator JENGSZ. 
Subroutine PERFM is used to compute take-off performance. Engine weights 
are computed by a call to subroutine ENGWGT (described in Volume V) . Subrou- 
tine TPALT (Section 1.1.3.15) is used to obtain atmospheric properties. Flap 
setting are determined by a call to subroutine APPFLP which is described in 
Section III. 1.4.4. Configuration drag is determined by a call to subroutine 
DRAG (Section III. 1.2. 2). Where engines are sired in turning flight, 
subroutine TURN is used (Volume VI) . 

A detailed flow chart for subroutine ENGSZ is presented in Figure iv.3.2. 


rv-3 

6 



BNGSZ 



JT0=0 

JRC=0 

JTR=0 

PAMB^PO 

TAMB=TO 

QA^^B* . 7*PAMB*EM*EM 

SMN=EM 

HN=H 

ASON=49.1*SQRT(TAMB) 

V=SMN*ASON 

wi;r=wg*rwcrto 

RiiUi:R=Ri;Li 
! Ql)=Q/'; B 
I IIRGOURG 


FIGURE IV. 3. 2 ~ DETAILED FLOWCHART , 
SUBROUTINE ENGSZ 
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ENGSZ 


2 
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ENGSZ 


3 


©- 


CDNAC-0. 

DRGNAC-0, 


c 


1 


KNAG 


c 




- iDRG-DRC+DRGNAC | 


JTO 


C 


> 



1002 


JRC 




C 


JTR 


> 


fcA?!" liNG I Nt: ( FNX , SFCX , SFNX , WFX , FAR . POX , TOx7eMx7hX 1 


R(:WACX=PCWAC 


C 


JTR 




->—{99 


K 

NAG 


< 

> 

< 

ICRU=1 

YDRG=DRG 

YCLB=WGR*RCCRU/V/60. 

i;npp=hnp 


' 



[call tINGINE(FN,SFC,SFNCR,WFR,FAR,PAMB,TAMB,SMN,HN,3) | 




IGRU=0 

pcwacy=pcwac 

WA= (DRG+WCR*RCCRU/ (V*60. ) )/ENP/SFNCR 
T2=TAMB* ( 1 . + . 2*SMN*SMN 
P2=PAMB* ( (T2/TAMB) **3 . 5) *PR 
WAC=WA*SQRT (T2/S 1 8 . 7) / (P2/21 1 6 . ) 
WASLS 1= WAG/ PAWACY 


> 

f 






JE 

NGSZ ^ 

> 

' 

1 WASI.S=WASLS1 

> 

f 


=4 


C wa ^,sx;;^3^ ^ 


WASLS 


> 


<WASLSX. 


4 wasTJs^waslsx 
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ENGSZ 


WRITE(6)WASLS, 

ERRT0,ERRM1 













ENGSZ 


C 


IRQMAX-5 n 

I 

y 


«i 


c 


c 


- |h-hport I 


IRQ 


I 


3 ^ 


- iH«HPORT+250.1 


C 


IRQ 


> 


H«HPORT+150g.) 


IRQ 


C 




-[h=hport3 


IRQ 


> 


|h=HPORT 




— < =3 

'' 

CALL TPALT(H,GALT,PO,FKALT,TO,G,XKV) | 

' 

f 


P0=PU*144. 

TO=TO+TDELTO 

RHO=PO/TO/53.32/G 



VSTO=SQRT(WGS/RHO/ GLMXTO/ .5) 
VROT=VSTO+10./.5929 i 

VLOF=VROT+5./.5929 ! 

VCLIMB=VLOF 
RCRQ-1. 

IDRAG=2 
ENf,OlJT=l. 


± 


t VCLMB»1.15*V"ST0^1 ^ 


VSTO=SQRT (WGS/RHO/CLMXTO/ . 5 ) 
VROT=VSTO+10./.5929 
VLOF=VROT+5./.5929 
VCLMB=VL05+45. 

I 


VCLMB 


^ <(1.15*VSTO) 


RCRQ=(VCLMB*60.)*.O24 

IDRAG=5 

ENGOUT-1. 
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bngsz 



VST0=SQRT(WGS/KH0/C 

VST0KT=VST0*.5929 

VCLMB*1.3*VST0 
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VSLD»SQRT(WGS/RHO/CLMXLD/ .5) 
VSLDKT«VSLD*.5929 
VCLMB=1.3*VSLD 
RCRQ=5.75*VSLDKT 















ENGSZ 



















IV-3 

17 














ENGSZ 


11 












ENGSZ 


12 
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ENGSZ 


14 











rv.3.1.3 Siibroutlna NACDG, Nac<ll« Subrcmtin>« Nl^CEXS cooaput*s 


lo«c«s during anglns sizing and parformanca calculations. Tha 
mathodology is discussed in Section IV. 1.1. 3. A detailed flowchart for 
subroutine NACDG is presented in Figure IV. 3. 3. No other subroutines are 
called by NACDG. 
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NACDG 



RETURN 


FIGURE IV. 3. 3 - DETAILED FLOWCHART SUBROUTINB NACDG 
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XV. 3. 2 Prop«ll«r Driven Sngin* Routin** 


IV, 3. 2,1 Subroutlna BNGINB, PropaillT Drlvam gnglna> R»rfonnanc« at 
3paclfl«d Flight Condition . Thi« routin* cctnput«« pi*ton ®ngin«, rotary 
combustion engine, or gas turbine engine pecforoance when matched to a 
specified propeller. Methodology is discussed in Section IV. 1,2. 2. Engine 
sizing options are exercised through the indicators KENG and KODE. The 
indicator NTYE controls engine type selection by calling either subroutine 
PWRPLT, (Section IV. 1.2. 3), for piston engine performance or subroutine 
TURBEG (Section IV. 1.2. 4), for turboprop engine performance. Atmospheric 
properties are obtained from subroutine TPALT (Section 1.1.3.15). Propeller 
characteristics are determined through subroutine ENGDAT (Section IV. 1.3). 

A variety of utility subroutines are called in the engine performance calcu- 
lations including BXV (Section I. 1.3. 4), MAXMHW, (Section 1.1.3.11), ITRMHW, 
(Section I. 1.3. 8), and MAXBND, (Section 1.1.3.10). 

A detailed flow chaurt for sxabroutine ENGINE is provided in Figture IV. 3.4. 
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EMGINE 



















ENGINE 


2 


XN-XNMAX*PCRPM*GR 

EFFP«1. 

ERRD-1. 

F*1.05 

FF-1.01 

JC«0 

JX=0 













ENGINE 


S 

































ENGINE 


5 
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BNGINH 


6 



RETURN 



CALL P WRP LT CHPM , 1 IPMS LS , HP WR , HP AVLB , PCPOWR , PCRPM , DE LTA , 


RTHLT , H , KSPCHG , BSFC , KODE) 







ENGINE 


KODE-8 

WF«BSFC*K }WR*HPMSLS 
TSFC-WF/THPROP 


■4 or 6 


PCRPM=1. 

F».9 

ERROR»999. 

JC»0 

JX=0 

KWRITY-0 


TMREQ=THRUST 

THPROP=THREQ-FN 

TPROPI=THPROP/(l.-FT) 


ERRMl=ERROR 


KWRITY-1 


CALL LNCDAT (PCRPM, GR , DPROP .TPROPI , SHP ,EFFDI , VKTS , 
RRHO , KODE , I ERROR , ENP , PO , TO , AF , BL) 


THPROP=TPROPI*(l.-PT) 

EFFP=EFFPI*(.1-FT) 

CT=CT*(1.-FT) 

HPRW=SHP 


CALL PWRPLT (llPM.llPMSLS , HPWR.HPAVLB, PCPOWR, PCRPM, 
DELTA , RTHET ,11 , KSPCHG , BSFC , KODE) 


EKROR= (S11P-11PAVLB)/I1PAVLB 


KWRITE 


WRITE(6)PCP0WR, 

PCRPM.SIIP.IIPAVLB 

ERROR 


ABS TERROR) 


<.001 












ENGINE 















ENGINE 


9 
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ENGINE 


10 









ENGINE 


11 











ENGINE 


12 
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ENGINE 


13 



RETURN 
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ENGINE 


14 





IV-3 

38 











IV. 3. 2, 2 Subroutln« gNGDAT, Prop»llar Qxttra.ctaLglgti.Gg . This routine 
controls propeller performance, cost, weight, and noise computations. 
Siibroutine PERFM (Section IV. 1.3) provides performance calculations! subrou- 
tine ZNOISE (Section IV. 1.3. 5), provides noise calculations; subroutine 
COST (Section IV. 1.3. 2) , provides costs; subroutine WAIT, (Section IV. 1.3.4), 
provides weight calculations. Gearbox characteristics are computed through 
subroutine GEARBX (Section IV. 1.3. 3). The utility routines BIV (Section 
1. 1.3. 4) and ITRM, (Section 1.1. 3. 7), are also ffinployed by subroutine ENGDAT. 

A detailed flow chart for ENGDAT is provided in Figure IV. 3. 5. 
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ENGDAT - PROPELLER VERSION 
















ENGDAT 


2 


TIPSD= . 05236*RPM*DROT 


<11 


XNOE=ENP 

FC=SQRT(518.7/T0) 


CALL ZN()ISE(BL.I)ROT,TIPSPD,VKTS,SHP,DIST,PNL,FC,XNOE) 


ZJI=5.309*VKTS/T1PSPD 
FC=SQRT(518.7/T0) 
ZMS(1)=.00152*VKTS*FC 
ZMS (2) =TIPSPD*FC/1120. 


CP= 1 0 . E 1 0*R0R0/ (2 . *TIPSPD**3*DR0T**2*6966 . 


CALL PliRFM ( KP1-;RFM , CP , ZJ I , AFT , BL . CLl , CT , ZMS , LIMIT) 


CALL P1;R1M(4, CP, Z.I1, AFT, BL.CLI.CT, ZMS, LIMIT) 


S1!P=CP*2*TIPSPI)**3*I)ROT**2/(10.E10*RORO)* 09966. 


THRUST=CT*TI PSPD**2*UR0T/ (1.51 5E06*R0R0) *364 . 76*XFT 

BLANG=BLLLL 

EFFP=CT*ZJI/CP 
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ENGDAT 


3 















ENGDAT 


4 



rv-3 

44 










IV.3.2.3 SuJaroutln* BWGSZ, propallaar Drivan ilgtog . ttibroutin® 

EKK3SZ controls th« propellar drivan angina aising calculations. Tha angina 
is sisad to jaaat cruise/ taka-off, esjid/or climb raquiraaiants using the 
aaathods of Section IV. 1.2.1. Subroutines called by ENGSZ include ENGINE 
for angina/propallar matching (Section rv.1.2.2)} PERFM for propeller perfor- 
mance (Section IV. 1.3.1); i^PFLP for flap setting (Section III. 1.4. 4); 

DRRG for configuration drag (Section III. 1.2. 2); TP2UuT for atmospheric 
characteristics (Section 1.1.3.15); ENGWGT for engine weights (Volume V); 
and the utility routines ITRMHW (Section I. 1.3. 8) . 

A detailed flow chart for subroutine ENGSZ is presented in Figure IV. 3. 6. 
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ENGSZ 


c 


NPC 


c 


KWRITE 


KFPTCH=0 


> 


> 


5<-l 


RETURN 


c 


I DC 


J - 


5^9 


WRITE (6) 


c 


JENGSZ 


> 


n,2 


XTORQ»99999. 


IANGLE»0 

XXDLRC«-99999. 

JNAC*0 

JTO*0 

JRC=0 

KCRV*0 

SMN=EM 

HN=H 

ASON*49.1*SQRT(TO) 

V=EM*ASON 

WCR=WG*RWCRTO 

POCRU^PO 

TOCRU=TO 

QCRU= . 7*P0CRU*EMCRU*EMCRU 

RI;LICR=RELI 

QO=QRCU 

DRG=URGC 


c 


NSC 


> 

f 

1 IFPTC 

o 

11 



i 


> 




'JFPTCH=0’1 

NTYP ') - 


c 


BLANG 


C 




JFPTCH-1 


JFPTCH 


3 


BLANG 






c 


JFPTCH 


> 


?^0 


KNAC 


> 


=2 


-5^ (18 


V HO 

T AWnt 1 

f * 



-© 


FIGURE IV. 3. 6 ~ DETAILED FLOWCHART, SUBROUTINE BNGSZ - PROPELLER 

DRIVEN ENGINE SIZING 
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ENGSZ 


2 


































ENGSZ 


4 
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ENGSZ 


6 



VST0-SQRt(WGS/RH0/CLMxW.5}^ 


/ > 

VR0T-VSTO+10./.5929 

VL0F-5./.5929 

VCLIMB=VLOF 



RCRQ«1. 

IDRAG=2 


1 

ENGOUT-1. 1 





63}— 


VSL11=SQRT(WGS/RH0/CLMX1)/ . 5) 
VClJ‘1B=1.5*VSLn 
RCRQ=(VCIWB*60.)*.021 
inRAG=8 

1 

: , 

CAM. 

l’l'Fl,P(Gi.MXU),I)FLAl', lOK} 


CZEZI)-^ 


Idrag»3”] 

i« momr 


r;Nc:ouT=i 
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ENGSZ 


7 


VSLD-SQRT (WGS/RHO/CLMXLD/ .5) ' 


VCLMB-1.3*VSLD 


RCRQ« (VCLMB*60.)*.032 


IDRAG»3 

> I’ 

ENGOUT=0. 



IRQM;\X=6 




H»HPORT 1 


( > 


-2 


c 


H=HPORT+250. 


IRQ 


3 


■3 


-> I H=5000. 


( > 


H=HPORT 


zr 


( > 


= 5 


T^s'ooor'i 


C ) 


= 6 


INIirORT 



• - - - 

< 1 

CAl.l. 

IPALI (11 ,GALT ,PO , I-KALT ,T0 ,G , XKV) 




r0=l’i.*144. 

'ro=TO+TlH;LTO 

RJIO=ini/TlV53.32/G 


■n 





VSTO=SQRT (WGS/RHO/GIiiXTO/ . 5) 


VST0=SQRT ( WGS/ RHO/CLMXTO/ . 5 ) 

VCIJ‘1B=1.2*VST0 


VSTOKT=VSTO*.5929 

RCRQ=50. 


VCLMB=1.2*VSTO 

II)IG\G=2 


RCRQ= . 046*VSTOKT*VST0KT 

r.NGOUT= 1 . 


IDRAG=5 

I 

ENGOirr= 1 . 


70l 
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ENGSZ 


CLMX-CLMXTO-DCLTO 
VS»SQRT(WGS/RHO/CLMX/l . 5) 
VSLD»SQRT(WGS/RHO/CLMXLD/ . 5) 
VSLDKT=VSLD*.5929 
VCLMB*1.2*VS 

RCRQ« ( . o79- . 1 06/ENP) *VSLDKT*VS 

IDRAG=1 

ENGOUT»l. 


VSLI)=SQRT(WGS/RflO/CLMXLD/.5) 

VSLI)KT=VSLD*.5929 

VCLMB=1.3*VSLD 

RCKQ= . 053*VSLDKT*VSLDKT 

II1R/\G=6 

i;ngout= I . 


CU4X=CLMXTO-I)CLTO 

VS=SQRT(KGS/RUO/CLMX/ .5) 

VSLD=SQRT(WGS/RHO/CLMXLD/.5) 

VSLI)KT=VSLn*.5929 

VCI.MB-:1.2*VS 

RCRQ=9.*VSLl)KT 

1DRAG=1 

r.NGOUT=0. 


VSLD=SQRT(WGS/RHO/CLMXLD/.5) 

VSLUKT=VSLI1*.5929 

VGLMB=1.4*VSLU 

RCRQ= . 092*VSLnKT*VSLDKT 

II)RAG=3 

i;ngout=o. 












ENGSZ 
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ENG5Z 


10 













ENGSZ 


11 
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ENGSZ 
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AN 


12 


riP~^TX«23 f 


WRITE (6)H1, VI, RCl, 
RCRQ1,CLR1,ZLQD1,H2, 

V2 , RC2 , RCRQ , CLR2 , ZLQD2 , 
H3,V3,RC3,RCRQ3,CLR3, 
ZLQD3 , H4 , V4 , RC4 , RCRQ4 , 
CLR4,ZLQD4 



WRITE (6) 
DFLAP 
















ENGSZ 


14 


HPMSLS 


<HPPAS 


(TN-TNPAS) 


WRITE (63 
TNjTNPAS 


-RETURN 


KWRITE 


WRITE(())RCRQ, 
liRRCL.ERRMl.DPROP, 
HPMSLS 


<.005 










ENGSZ 


15 


-1 
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ENGSZ 


16 
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IV. 3. 2, 4 gubroutijv CSaARBX/ <g#Rrb03£ Weight/ Qoxt, and Molsa . Subroutine 
curries out gssrbosc weight/ noise, and cost computations using the 
toethod of Section IV. 1.3. 3. No ottier subroutines are called by this routine. 
A detailed flow chart for subroutine GEAEIBX is presented in Figure IV. 3. 7. 
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FIGURE IV. 3. 7 - DETAILED FLOWCHART, 


SUBROUTINE GEARBX - GEARH O X WEIGHT , 
COST AND NOISE 
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IV. 3. 2. 5 gubroutlna PBRffM, Propallor Parfamanctt . SubroutJja* FEKTO 


oomputas propallar parformanca by tha jnathod of Saction IV. 1.: 
latlon* mainly involve interpolation in atorad data ixalng tha 
routine* BIV (Saction I. 1.3. 4) and UNINT (Saction 1.1.3.17). 
flow chart for aubroutina PERFM ia praaantad in Figura IV. 3. 8 


1 . 1 . Calcu- 
utility 
A datailad 


rv~3 



PERPM 



CALL UNINT(6,AFVALC1) ,AFCPC(1,K) ,AFT,AFCP(K) .LIMIT) 
CALL UNINT(6,AFVAL(1) ,AFCTC(1,K) , AFT, AFCT(K) , LIMIT) 



135j FIGURE IV. 3. 8 - DETAILED 

FLOWCHART, SUBROUTINE PERFM 
PROPELLER PERFORMANCE 

















PBRFM 


2 



DO 500 IBB=1,NBB 

I 


DO 300 K»NBEG,NEND 



DO 213 KL=NCLT,NCLTT 



CALL UNINT(NCLX(NNCLT) .CPCLI (1 ,NNCLT) ,XPCLI (1,NNCLT) , 
CPE1,PXCLI(KL) , LIMIT) 
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PERFM 


3 




c 


CPEl 




.049 


C > 


=5 


WRITE (6) CPEl, 
PXCLI(KL),L 


1 — ^ — ; < » 

( CPEl > 

i M - - 

WRITE (6)CPE1, 
PXCLI(KL),L 



c 


CPE I 




0915 


WRITE (6)CPE1, 
PXCLI(KL),L 


NLRPT=1 

NNCLT=NNCLT+1 


!15 CONTINUE 


c 


nc:e 


> 




|CA1 1, UN INT (4 , CCLI (NCLT) , P}^LI (NCLT) ,CLI , PCLI , LIMITrj 



jT{Tr=PX(:U (NC I 

^ 

lcp;:=cPE*PCLi 


CAI.l. UNINT{INN(K),CPATIG(1,K,L3,BLDANG(1,K3,CPE,BLL(K),L1MITJ 
CAU. UNINT(INN(K).BLDANG(1,K) ,CTANG, Cl , K, L) ,BLL(K) ,CTT(K) , I.IMIT) 


C 


LIMIT 


> 


=0 


-> (300, 



|nERPT=2 I 




|NNCLt=NCI?fl 


DO 250 KL-NCLT.NCLTT 
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' 

f 

CTA(1] 

CTAC21 

*CT 

=1.5*CT 












PERFM 


4 


DO 2600 KJ«1,5 


NFTX-KJ 

CTE1=CTA(KJ)*AFCT(K) 


' 

L 

|CALL UNINT(14,CTEC(1) ,BTDCR(1,L) ,CTE1 ,TBL, LIMIT) | 



CTr.i=c 

rLl+TBL*TFCLl (K) 


CALL UNINT(NCLX(NNCLT) ,CTCL(1,NNCLT) ,XTCLI (1 ,NNCLT) , 
CTL1,TXCLT(KI.) .LIMIT) 


nl;rpt=3 


C 


LIMIT 





591 


^ ZJJ(K) J - 


ZMCRT«ZMCRO(NNCLT) 

DMN=ZMS(2)-ZMCRT 


}cAL^ UN 1 NT (Tl ,ZJCL ( 1 ) , ZMCRL ( 1 , NNCI.tT, ZJJ ZMCRT, LIMIT) | 


DMN=ZMS(1)-ZMCRT 


lxn-T(Ki.) = i.ti 


-c 


I 


DMN 


3 



ri;i*TXCLI(KL)/TFCLI(K) 



1 CALL B I gilAU ( ZMMC , 1 , DMN , CTH2 , XFFT (KL) , LIMIT) j 

> 


}t:iAl lKJ) --CT-CTA{KJ)*XM TlKi.) j 


CTAUK.l) ^ 


KJ 


> 




( ^PS (CTA I ( KJ - 1 ) - CTA 1 ( KJ ) /C^ 

'' 


,001 


IV-3 

68 

















PERFM 


5 



CALL U.\ l NT 1 4 . Z. I J (NBlii I ) . BLL (NBEG) , ZJ I (NBLCO , 
BI,I,1,( IBB) .LIMIT 
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PERFH 


6 



DO 39b KL=NCLT,NCLTT 



































PERFM 


S 


DO 290 IL«1,5 


t 


CP-CPG(IL) 
CPE*CPG(IL) *AFCPE 





jCALL UNINT(14,CPEC(1),BLDCR(1,L), CPE, PBL, LIMIT)! 



CPL1=CPE*PBL*PFCLII 

NNCI.T=NCLT 



t 

DO 280 KL*NCLT,NCLTT 


|CALL UN INT (NCLX (NNCDn , CPClYc 1 ,NNCLT) , XPCLI ( 1 ,NNCLT) , CPEl , PXCLI (KL) , LIMIT)*) 















PERFM 


9 















IV.3.2.6 Subroutlna PMOYSf Drlvan Aircraft Nolaa Controlling 

Routlna . Svibroutina PNOYS controls the propeller driven aircraft noise 
calculations as discussed in Section IV. 1.3. 6. Routines called by PNOYS 
incliide subroutine ENGINE for engine performance (Section IV. 1.2. 2) f 
stibroutine GEARBX for gear box characteristics (Section IV. 1.3. 3); Subroutine 
ZNENG for engine ncjise characteristics (Section IV. 1.3. 7) j subroutine TPALT 
for atmospheric properties (Section 1.1.3.15), and subroutine ASPEED. 

A detailed flow chart for subroutine PNOYS is presented in Figure IV. 3. 9, 
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PNOYS 



FIGURE IV. 3. 9- DETAILED FLOWCHART, SUBROUTINE PNOYS 
PROPELLER DRIVEN AIRCRAFT NOISE CALCULATIONS 
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PNOYS 


2 




Rli'llIRN 


rv~3 

76 










IV. 3. 2.7 Subroutlnai gIfSPLg, flaton Bn-glna "Bomr and Fual Flow . 
Svibroutina PWRPLT computes piston angina power and fuel flow by tha mathod 
of Section IV. 1.2. 3. Tha only svibroutina called by PWRPLT is tha utility 
routine ITRm. A detailed flow chart for subroutine PWRPLT is presented in 
Figure IV. 3. 10. 
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PtfRPLT 


[call ITRLN (APCN, APCPWR , PCRPM, PCPWKM, 6) 


SU:MA=DHLTA/ ( (RTHLT) ** 2 ) 


=2. 5. 6J, 8. 9, 10,11 


■1,3,4 


I Hl’MiLS=HPWR/PCPOWR 


*5 or 6 or 9. 


PCPOWR»HPWR/HPMSLS 


»2 or 7 


HPWR=HPMSLS*PCPOWR 


KSPCHG 


iCAI.I. ITRLN(A1I.AW1P.H,RUPC0R,6) 


llPMi:OR=llPMSLS*WlPCOR 
1 1PM= { IPMCOR* UP LTA* RTl H;T 



SIGMA 


HPM=1IPMSLS 


<SIGCRT 


1 1PM= ( (S I GMA- . 1 1 7 ) / (S IGCRT- .117)) *HPMSLS 


HPAVLB=11PM*PCPWRM 





<HPWR 



(HPWR-HPAVLB) /HPAVLB 


HPAVI.B=HPWR 


HPWR=UPAVLB 

PGPOWR=HPWR/HPMSLS 


=11PAV1.B 


FIGURE IV. 3.10 -- DETAILED FLOWCHART 
SUBROUTINE PWRPLT, PISTON ENGINE 
POWER AND FUEL FLOW 


HPAVLB=IiPWR 

iipm=hpavi.b/p{:pwrm 










PWRPLT 


2 



( KSPCHG ~ ) 


J^l 


c 


HPMCOR-HPM/DELTA/RTHET 

HPMSLS»HPMCOR/RHPCOR 


SIGMA 


> 


<SIGCRT 


HPMSLS=HPM/( (SIGMA-. 117) 
/(SIGCRT-.117)) 


I UPMSl,S=llPM I 

1 


I PCPOWR=HPWR/HPMSLS 1 

1 

' 


■fniwawi 

i|2y|puV 




^ DISPL»2.*HPMSLS1 


P ISPL=792000 . *HPMSIS/BMEP/XNM\X j 


>.78 

>|TP ... 






GPH=(.0072+ 

.0365*PCPOWR)*DIS1’[,+O.I8 | 

1 


=8 

iWii 



' 


r 


C 


PGPOWR 


C 




f GPH» (7 o0338+ . 




0355*PCPOWR)*DISPl 


PCPOWR 


C 




78 


'spi,7o.,y*| 


GPtlN((Tp8^^CP0WR-rbF 


PCPtIWR 




= .78 


>t^ 


Hsic-i.. i)*i:i’ii/ih:i’owr/iipmsls 


C 


xnmax 


3 


•2700. 


->• RETURN 




j PHOB=(. 000045) *XNM/\X)-. 01447 
1 BSI-t;=BSl'C* ( 1 . +FHOBH) / 1 . 1 1 

> 

i IV- 3 


Rl'.TURN 
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1V,3,2«8 jEg^oatlna TOKBgQy y'^bofacop Bogin# P«gga t 3» yt i>^ » Saih^coiitiB,* 


TOKBEQ coc^-tes the performanc© of turboprop engines by the wetbod of 
Bectioa IV. 1.2. 4. The only svibroutinas called are the utility routines 
BXV (Section I. 1.3. 4); ITRLN (Section I. 1.3. 7); and ITBMHW (Section I. 1,3,8), 
A detailed flow chart of TURBEG is presented in Figure IV. 3. 11. 
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TURBEG 



FIGURE IV . 3.11 - DETAILED FLOWCHART, SUBROUTINE TURBEG 
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TURBEG 


2 


CALL ITRLNCAXNCR,ATCR,XNCR,T4QT2M,8) 
CALL ITRLN (AXNCR, ATCR , XNCREF , T4T2RF , 8) 





















TURBEG 


3 


CALL B I V (HPCT5 0 , XNCR , T4QT2 , AXNCR , AT4QT2 , AHPC5 0 , 8 , 6 , NER) 

CALL BI V (HPC50R , XNCREF , T4T2RF , T4T2RF , AXNCR , AT4QT2 , AHPC5 0 , 8 , 6 , NER) 



ITAB=4 


HPCTAB=HPCT25+ CHPCT50-HPCT25) * (SMN- . 25) / . 25 
HPCREF=HPC25R+ (HPC50R-HPC25R) * (SMN- . 25) / . 25 


CALL ITRLN (AXNCR, ATRT) , XNCREF, T4T2SL, 8) 

CALL BIV(HPSLRF, XNCREF, T4T2SL, AXNCR, AT4QT2,AHPC0, 8, 6, NER) 


RHP=HPCSZ/HPCTAB 

HPAVLB=HPWR 

HPCMAX=HPCREF*RHP 

HPM=HPCMAX*DELTA2*RTHET2 

HPMSLS=HPSLRF*RHP 

.IJ=1 


RHP=HPMSLS/HPSLRF 

HPAVLB=HPCTAB*RHP*DELTA2*RTHET2 

HPCMAX=HPCREF*RHP 

HPM= 1 II’CMAX * DE LTA2 * RTHET2 


ERRMUERRHP 


>PCNCMX 


ERR1IP= (HPAVLB-HPWR)/HPWR 












TURBEG 


4 
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TURBEG 


5 


















TURBEG 


6 

















TURBEG 


7 


WRITE (6) 


WRITE(6)ERRHP, 

1 

WRITE (6) JC, 

ITAB.SMN, 


ERRM1,JC,JC1,HPM, 

1 

JCi,T4QT2, 

XNCR.T4QT2 


HPWR/M 

1 



c 


KWRITE 


3 ^ 


or *9 


WRITE ( 6) RHP , SMN , T2 , RTHET2 , 
HPAVLB , XN , HPSLRF , XNCR , XNCREF , 
PCN , PCNCMX , PCNCR , T4QT2M , T4QT2 , 
T4T2RF,HPM,HPMSLS 


RETURN 
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IV. 3. 2. 9 Btibroutina Propallar Walght . ixJarcmtina WMT ocxciputa* 

jwropallar vraighta by tha mathod of Saction IV, 1.3. 4. Tha indicataar rWTCON 
datarminaa whici' ona of flva aata of aquationa ara uaad to predict 1970 
and 1980 propallar weights. No other subroutines ara called by WAIT. A 
detailed flow chart for aubroutina WAIT is provided in Figrura IV, 3. 12. 
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WAIT 
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WAIT - PROPELLER WEIGHTS 









i:V.3.2.10 SuJaroutirii^ gHBNg^ Saqlnd Noi«» . Tlid« routine ooiaput.** piston, 
rotairy, «nd turboshftft •ngin# noiss ch*ract«ri»tic« by th« meithod of Soction 
rv. 1.3.7. En»gin« typa is salactsd by th« indicutor ITYPK. Tha only aub- 
routina oallad by ZNENG la tha utility routina UNINT (Saction 1.1.3.17). 

A datailad flow chaxt for avibroutina ZNENG ia praaantad in Figura IV. 3. 13. 
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ZKENG 


XNOE-NOE 

X-ALOGIO(SHP) 

XX-10.*ALOG10(XNOE) 


ITYPE 


X=10.*X+XX 

F.Nni»28.5+X 

EBAI«16.3+X 

ENDC«51.6+X 

EBAC-38.8+X 

XNCYL=4. 


- 210 . 


XNCYL-6. 


XNCYL-8. 


H2FF=XNMAX*PCRPM* . 50*XNCYL/b0. 


CALL l)NINTt4,IlX:i'H{n ,XXUn .liZFF.A.LIMIT) 
CALL IININT(4,HXZFF(1) ,XX2(1) .MZFF.B, LIMIT) 


LIMIT 


= 0 

- — - >-• 4 


WRITE(b) 
HZFF _ 


ENm;=A+x 

i;bae=b+x 


X=10. 

*X+ 

XX 

ENPI = 

47. 

3+X 

e;bai = 

3b. 

7+X 

ENPC= 

4P 

+X 

E.BAC= 

28. 

1+X 

ENPE= 

72. 

9+X 

E.BAE> 

b2. 

2+X 


X-7.5*X+XX 
EPNDB-67.1+X 
EDBA=54.7+X 
F.N1U = 0. 

F.BA1 = 0. 
ENDC-O. 
EBAOO. 

IiNPE = 0. 
F.BAE=0. 


EPNDB«10.*ALOC.10(10.**(.1*ENDI) + 10.**(. I*F.NDC) + 10. ** ( . 1*ENDE) ) 
EnBA»10.*AL(H:i0(10.**{. 1*EBAI)+10.**(. 1*EBAC)+10.**{. 1*EBAE)) 


KK rn; ( b ) L NP I , i; BA I , LNl )C . KOH , 
LBAC . F.NDli , EBAE . liPNPB . EPBA 


Fi r.URK TV. 3. 13 - DPTAILED FLOWCHART 
SUBROUTINE ZNENG - ENGINE NOISE 







IV. 3. 2. 11 Subroutine ZNOIgB^ ggo^llar Noise . This routine cxjmputa* 
prqE>aller generated noise by the method of Section XV. 1.3. 5. The only 
subroutine called is the utility routine BILINE (Section I. 1.3.1). A 
detailed flow chart for subroutine ZNOISE is presented in Pigura IV. 3. 14. 
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ZNOISE 



WRITH(b)XNOE,VKTAS, 

mST,XLl,XL2.XL3,XL4, 

PNI.D,SPL,SPLX 


RETURN 


FIGURE IV. 3. 14 


DBTAILBD FLOWCHART, SUBROUTINE 
ZNOISE ~ PROPELLER NOISE 



















IV. 3. 2, 12 Svibroutliva COBV, Propallar Costa . This routiaa oomiputaB 
propallsr ooata by tha mathod of Saction IV. 1.3. 2. Both If/O or 1980 coat 
aatiaiataa may ba mada. No aubroutinas ara callad by coat, A datailed flow 
chart for atabroutina CXDST la praaantad in Figura IV. 3, IS. 
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COST 


ICON«WTCON-«-.01 I 


lENT CK70 


nniSBi 

+ 

> 

r 

' 

miiiiiiiiiii 

CCLF1 = 
CCLF= 

=3.2178 

1.02 

CCLP1«CLF1 

CCLF=CLF 


r ' 

' 


RETURN RETURN 


1 


CCK70- 



=ZFFAC ( 1 , ICON) * (3 . 0* BLADT** . 75+ZEFAC (ICON) ) 

1 ,M||— I.I.... 1 1 >.I.I II ■■—.WU - --..1 I-.L.IU ■ 

1 

^ 4. 


CK80 } f j CCK80=CK80 

>— 

• 

^ 0 

1 CCK80 

=ZFFAC{2 , ICON) * (3. 0*BLADT* * . 75+ZEFAC (ICON) ) | 


t 



RETURN 


FIGUKB IV. 3. 15 - DETAILED FLOWCHART, SUBROUTINE CX^ST— 

PROPELLER COSTS 
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IV.3.2.13 Subroutiniaa BMGDTT, BNGDTl to ENGPT7 

Thaaa routinaa provica propulaion angina charactariatica for variotia 
tairbojat anginaa. A datallad flowchart ia prasented for ENGDTl. Tha 
anginea daacribad by thaae routinea eo:a listed in Section IV. 3. 1.2. 
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BNGD.KX*- TURBOJCT VERSION 


THETA2-T2/S18.7 

RTHET2-SQRT(THETA2) 

DELTA2*P2/2116. 


> 

U!i_ 

KENG ^ 

t 

1 PCNCMX=100./RTHET2 


c 


PCNCMX 


C 


> 


>105 


-> 1pCNPCMX»104.98 1 


KHNG 




,5 


PCNCR=PCNCMX 


c 

c 




KPNG 


3 -^ 


.4 


PCNCR= PCNCMX *FTHROT 

""rr'Z 


KENG 


> 


=3.6.7 


PCNCR=98.0 


ITAB=1 


CALL B I V f FNCR . SMN . PCNCR . ASMN . APCNCR . AFNCR .5.8, NER1 


C 


NER 


> 




WRITE f6)ITAB. 
SMN, PCNCR 


RETURN 


ITAB=2 


CALL B I V ( WFCR , SMN , PCNCR , ASMN , APCNCR , AWFCR ,5,8, NER) 
CALL I TR LN (APCNCR , AWACR , PCNCR , WACR , 8 ) 


PCWAC=WACR=WACR/43.9 

FN=FNCR*DELTA2 

WF=WFCR*nELTA2*RTHET2 

WA=WACR*I)ELTA2/RTHET2 

SFN=FN/WA 

SFC=WI7FN 

FAR=WF/WA/3600. 


I 


WF=600,*DELTA2*RTHET 

WA=13.9*DELTA2/RTHEt 

FN=195.*DELTA2 

PCWAC=13.9/43.9 

FAR =WF/WA/ 3600. 

SFN=FN/WA 


— -f— 
RETURN 


FIGURE IV. 3. 16- DETAILED FLOWCHART, SUBROUTINE ENGCT3. 
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